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experimental  flow  rate,  at  a given  pressure  drop}). 

Empirical  models  have  also  been  developed  for  the  NSVR  configuration 
with  flow  in  the  low  resistance  direction,  and  for  the  SVR  configuration.  Thes 
models  have  been  developed  by  using  only  a limited  amount  of  data  (8  geometries 
Their  usefulness  is  correspondingly  limited. 

Approximate  analytical  models  have  been  developed  for  the  symmetric 
and  non-symmetric  vortex  resistors.  Two  steady-state  analytical  models,  one 
incorporating  a one-dimensional  and  the  other  a two-dimensional  approximation 
to  the  flow  field  in  the  vortex  chamber,  have  been  developed  for  the  non- 
symmetric  vortex  resistor.  A one -dimensional  model  has  been  developed  for  the 
symmetric  vortex  resistor  configuration.  These  models  exhibit  a maximum  error 
of  15  percent  of  reading  over  the  range  of  the  experimental  data  measured. 

The  feasibility  of  designing  temperature-insensitive  and  compensat- 
ing vortex  resistors  in  practical  sizes  has  been  established.  (Compensating 
vortex  resistors  cancel  the  temperature  sensitivity  of  conventional  fluid 
resistors . ) 


Temperature  sensitivities  ranging  from  +1.17x10 
-2.02x10"^  m^/s-  C have  been  measured  with  the  a 
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0 m /s-C  (almost 
available  vortex 


zero)  to  -2.02x10"'  mJ/s-  C have  been  measured  with  the  available  vortex 
resistor  configurations.  The  temperature  sensitivity  of  a comparable  short- 
tube  orifice  is  of  the  order  of  +1.47x10“^  m^/s-  C. 

A linearized  dynamic  model  based  on  a one-dimensional  flow  approxi- 
mation has  been  developed  which  incorporates  -mplicitly  a pure  time  delay  and 
a first  order  lag.  The  model  can  be  used  to  predict  the  response  of  the 
resistor  output  pressure  to  step  (with  finite  rise  time)  changes  in  the  input 
pressure.  Analytical  predictions  are  in  fair  agreement  with  the  experimentallj 
observed  dynamic  responses. 

Transient  responses  to  a terminated -ramp  input  measured  on  a typical 
geometry  indicated  that  the  rise  time  is  13  milli-seconds.  The  predicted 
rise  time  using  the  analytical  model  is  9 milli-seconds. 

Noise  measurements  made  on  three  typical  NSVR  geometries  (6.35mm 
vortex  chamber  diameter)  indicate  that  the  signal-to-noise  ratio  (for  signal 
levels  up  to  1722  kN/m^  pressure  drop)  of  a vortex  resistor  operating  in  the 
laminar  flow  regime  is  better  than  900. 
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NuMENC  LATURE 

o 

cross-sectional  area  of  the  loading  orifice,  (m  ) 
width  of  the  inlet  orifice,  (m) 

discharge  coefficient  of  inlet  orifice,  (dimensionless) 
discharge  coefficient  of  the  central  orifice,  (dimensionless) 
discharge  coefficient  of  the  exit  orifice  of  SVR,  (dimensionless'! 
discharge  coefficient  of  the  loading  resistor,  (dimensionless) 
height  of  the  vortex  chamber,  (in) 
length  of  the  capillary  tube,  (m) 
length  of  the  inlet  orifice,  (m) 

length  of  the  central  orifice,  (m)  2 

static  pressure  at  the  upstream  of  the  vortex  resistor,  (kN/m  ) 
static  pressure  at  the  downstream  of  the  vortex  resistor,  (kN/m2) 
static  pressure  at  the  upstream  of  the  central  orifice  [see  equation 
(B.13)],  (kN/m2) 

static  pressure  at  the  outer  periphery  of  the  vortex  chamber,  (kN/m2) 
static  pressure  at  the  exit  of  the  vortex  chamber,  (kN/m  ) 2 

static  pressure  at  the  downstream  of  the  central  orifice,  (kN/m  ) 
static  pressure  at  the  inlet  of  the  outlet  orifice  of  the  SVR,  (kN/m^) 
volumetric  flow  rate,  (m^/s) 
resistance  of  fluid  resistor,  (kN-s/ni  ) 
radial  coordinate,  (m) 

Reynolds  number  based  on  the  diameter  d,  (dimensionless) 

Reynolds  number  based  on  a diameter  dg,  (dimensionless) 

radius  of  the  vortex  chamber,  (m) 

radius  of  the  central  orifice,  (m) 

temperature  of  the  operating  fluid,  ( C) 

time  delay,  (s) 

first  order  lag,  (s) 

radial  velocity,  (m/s) 

radial  velocity  at  the  inlet  of  the  vortex  chamber,  (m/s) 
radial  velocity  at  t=T2  in  the  vortex  chamber,  (m/s) 
radial  velocity  at  r=r£  in  the  diffuser  chamber,  (m/s) 
radial  velocity  functions  [see  equation  (D.7)] 
tangential  velocity,  (m/s) 

tangential  velocity  at  the  inlet  of  vortex  chamber,  (m/s) 
tangential  velocity  at  the  inlet  of  the  central  orifice  (m/s) 
tangential  velocity  at  the  downstream  of  the  central  orifice, 

(m/s) 
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characteristic  length  used  for  shear  stress  calculations  (see 
Appendix  B) , (m) 
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Appendix  B),  (m) 
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NOMENCLATURE  (cont'd) 

incremental  change  in  the  variable  in  the  parenthesis  (eg. 
pressure  drop  across  a resistor,  (kN/m^)  1 

change  in  volumetric  flow  rate,  (mJ/s) 
boundary  layer  thickness , (m)  2 

a dimensionless  parameter  defined  by  24p.  r^/h  I Uj| 
absolute  viscosity  of  the  fluid,  (kN-s/m^)  ' ' 

kinematic  viscosity  of  the  fluid,  (m^/s) 
fluid  density,  (kN-s^/nA)  2 
tangential  shear  stress,  (kN/m  ) 

tangential  sheajj  stress  on  the  annular  wall  of  the  vortex 
chamber  (kN/ni  ) 

tangential  shear  stress  in  the  vortex  chamber,  (kN/nr^ 
tangential  shear  stress  in  the  central  orifice,  (kN/ni  ) 


1.  INTRODUCTION 

1.1  Background 

In  fluid  power  and  fluidic  systems  resistors  are  used  for  var- 
ious applications  such  as  impedance  matching,  flow  rate  measurement  and 
control.  Orifices  and  capillary  tubes  are  the  most  common  fluid  resist- 
ors. Figures  1 and  2 show  typical  steady-state  characteristics  of  a 
capillary  tube  and  a short-tube  orifice  respectively.  Hydraulic  fluid 


Figure  1.  Steady-state  characteristics  of  a capillary  tube 
(assuming  fully-developed  laminar  flow). 


Figure  2.  Steady-state  characteristics  of  a short-tube 
orifice  with  sharp-edged  entry. 
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is  assumed  to  be  the  operating  medium  here  and  throughout  this  report. 
Normally  AP  is  considered  the  independent  variable  (plotted  on  abscissa) 
and  Q the  dependent  variable  (plotted  on  ordinate).  Resistance  is  de- 
fined as  the  inverse  slope  of  the  Q vs.  Ap  curve  at  an  operating  point, 
i.e.,  Rs^(AP)/5Q.  Since  most  fluid  resistors  are  nonlinear  (R  depend 
on  the  operating  point),  it  is  convenient  to  describe  the  steady -state 
behavior  in  terms  of  the  Q vs.  AP  relationship  or  "characteristic  curve" 
instead  of  the  "resistance."  Resistor  geometrical  parameters,  fluid 
thermodynamic  properties  (density  and  temperature  or  viscosity)  determine 
the  steady-state  characteristics.  The  term  "temperature  sensitivity," 
used  in  this  report,  describes  the  variation  of  the  steady-state  flow 
rate  through  the  resistor  due  to  fluid  temperature  variation,  i.e., 
evaluated  at  an  operating  point.  For  fluid  temperature  ranges 
normally  encountered  in  practice,  the  temperature  sensitivity  of  a 
resistor  can  be  substantial. 

In  many  applications  the  variation  of  flow  rate  (at  a constant 
pressure  drop)  with  the  temperature  is  undesirable.  For  example,  if  the 
pressure  at  the  supply  port  of  a hydraulic  proportional  amplifier  is  held 
constant  (normal  operation),  the  supply  flow  rate  may  vary  substantially 
due  to  the  sensitivity  of  the  orifice  resistance  to  temperature  varia- 
tions. This  variation  affects  the  gain  characteristics  of  the  amplifier 
and  thereby  alters  steady-state  and  dynamic  behavior  of  the  associated 
system.  In  this  case  it  is  desirable  to  use  a compensating  resistor 
to  minimize  the  temperature  sensitivity  of  the  supply  port. 

Vortex  resistors  can  be  designed  to  have  a positive,  zero,  or 
negative  temperature  sensitivity.  Capillary  tubes  and  short-tube  orifices 
always  exhibit  positive  temperature  sensitivities.  Properly  designed 
vortex  resistors  may  be  used  with  capillary  tube  and  orifice  type  resis- 
tors to  achieve  temperature  compensation. 

Figures  3 and  4 show  schematics  of  two  types  of  vortex  resis- 
tors: the  Non-Symmetric  Vortex  Resistor,  or  "NSVR"  (also  called  the 

vortex  diode),  and  the  Symmetric  Vortex  Resistor,  or  "SVR."  Analysis  of 
the  fluid  mechanics  of  the  vortex  resistors  reveals  that  the  resistance 
elements  (orifices,  vortex  chamber  and  diffuser  chamber)  can  be  divided 
into  two  groups.  One  group  (vortex  chamber  and  central  orifice)  contri- 
butes a negative  temperature  sensitivity.  In  contrast,  the  second  group 
contributes  a positive  temperature  sensitivity  effect  and  tends  to  cancel 
the  effect  of  the  first  group.  This  second  group  consists  of  the  inlet 
orifice  and  the  sudden  expansion  diffuser  in  the  case  of  the  NSVR.  In 
the  SVR  configuration,  the  inlet  and  exit  orifices  together  with  the 
diffuser  chamber  form  the  positive  temperature  sensitivity  group.  Figures 
5 and  6 show  measured  steady-state  characteristics  for  two  NSVR  geometries, 
one  which  exhibits  a negative  and  one  a small  positive  temperature  sen- 
sitivity. When  the  positive  and  negative  effects  of  the  two  resistance 
groups  exactly  match  each  other,  the  resultant  vortex  resistor  is 
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temperature- insensitive.  Figure  7 shows  a case  in  which  the  temperature 
sensitivity  is  near  zero,  or  negligible. 
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Figure  7.  Measured  steady-state  characteristics  of  a NSVR 
with  negligible  temperature  sensitivity. 

For  a given  pressure  drop  and  a comparable  minimum  dimension 
(i.e.,  vortex  chamber  height  compared  to  orifice  width  or  diameter), 
the  vortex  chamber  can  offer  substantially  higher  resistance*  than  the 
simple  orifice  alone.  Figure  8 shows  the  pressure  drop-flow  rate  char- 
acteristics of  two  different  vortex  resistors  and  for  a simple  short- 
tube  orifice  which  has  the  same  geometry  as  the  central  orifice  in  the 
vortex  resistors.  The  increase  in  resistance  due  to  the  inclusion  of 
the  vortex  chamber  is  apparent.  Thus,  a given  resistance  can  be  achieved 
with  a larger  minimum  dimension  in  a vortex  resistor  than  in  a simple 
short-tube  orifice.  Sensitivity  to  contamination  due  to  solid  particles 
in  the  hydraulic  fluid  is  correspondingly  less  for  the  vortex  resistor. 


1.2  Scope  of  Study 

The  principal  objectives  of  this  study  were: 

1.  To  develop  and  validate  mathematical  models  which  describe 
*R  s h (AP)  /*Q 
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Figure  8.  Typical  predicted  steady-state  characteristics 
of  non- symmetric  vortex  resistors  and  a comparable 
short-tube  orifice  for  constant  temperature  operation. 

the  steady-state  and  dynamic  behavior  of  the  two  basic  types  of  vortex 
resistors,  i.e.,  Q = Q(AP,T,t). 

2.  To  determine  by  experiment  on  typical  geometries  the  noise 
characteristics  of  vortex  resistors. 

A secondary  objective  was  to  demonstrate  the  feasibility  of 
designing  temperature-insensitive  vortex  resistors,  and  temperature- 
sensitive  vortex  resistors  which  can  compensate  for  the  temperature  sensi- 
tivity of  a conventional  fluid  resistor. 

The  study  was  limited  to  the  following  specific  tasks: 

1.  Development  of  steady-state  empirical  models  for  the  two 
basic  types  of  vortex  resistors. 

2.  Development  of  approximate  steady-state  analytical  model 
for  the  two  types  of  vortex  resistors. 
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3.  Development  of  an  approximate  dynamic  model  for  the  non- 
symmetric  vortex  resistor. 

4.  Measurement  of  noise  characteristics  of  three  typical 
vortex  resistors  (non- symmetric  configuration). 

5.  Determination  of  maximum  and  minimum  steady-state  tempera- 
ture sensitivities  for  a set  of  geometrical  configurations  which  cover 
the  expected  practical  range. 


a 


1.3  Summary  of  Results 

An  empirical  model  has  been  developed  for  the  NSVR  configuration 
(for  flow  in  the  high  resistance  direction).  Experimental  data  taken 
on  24  different  geometries  and  at  temperatures  of  35  C to  57  C (with 
MIL-5606A  hydraulic  oil)  have  been  used  in  the  model  development.  The 
model  predictions  exhibit  a maximum  error*  of  30  percent  over  the  range 
of  the  experimental  data. 

Empirical  models  have  also  been  developed  for  the  NSVR  config- 
urations with  flow  in  the  low  resistance  direction,  and  for  the  SVR 
configuration.  These  models  have  been  developed  by  using  only  a limited 
amount  of  data  (8  geometries) . Their  usefulness  is  correspondingly 
limited. 


Approximate  analytical  models  have  been  developed  for  the 
symmetric  and  non-symmetric  vortex  resistors.  Two  steady-state  analytical 
models,  one  incorporating  a one-dimensional  and  the  other  a two-dimen- 
sional approximation  to  the  flow  field  in  the  vortex  chamber,  have  been 
developed  for  the  non-symmetric  vortex  resistor.  A one -dimensional  model 
has  been  developed  for  the  symmetric  vortex  resistor  configuration. 

These  models  exhibit  a maximum  error  of  15  percent  of  reading  over  the 
range  of  the  experimental  data. 

The  feasibility  of  designing  temperature- insensitive  and  com- 
pensating vortex  resistors  in  practical  sizes  has  been  established. 
(Compensating  vortex  resistors  cancel  the  temperature  sensitivity  of 
conventional  fluid  resistors.)  Temperature  sensitivities  ranging  from 
+1. 17x10'” nr /s-*C  (almost  zero)  to  -2.02x10"' nr /s-  C have  been  measured 
with  the  available  vortex  resistor  configurations.  The  temperature 
sensitivity  of  a comparable  short-tube  orifice  is  of  the  order  of 
+1.47x10"'  m ^s-*C. 


Q, 


*Error  = 100  x (experimental  Q - predicted  Q) /Experimental 
at  a given  AP. 
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A linearized  dynamic  model  based  on  a one -dimensional  flow 
approximation  has  been  developed  which  incorporates  implicitly  a pure 
time  delay  and  a first  order  lag.  The  model  can  be  used  to  predict  the 
response  of  the  resistor  output  pressure  to  step  (with  finite  rise  time) 
changes  in  the  input  pressure.  Analytical  predictions  are  in  fair  agree- 
ment with  the  experimentally  observed  dynamic  responses.  Transient 
responses  to  a terminated-ramp  input  measured  on  a typical  geometry 
indicated  that  the  rise  time  is  13  milli-seconds . The  predicted  rise 
time  using  the  analytical  model  is  9 milli-seconds. 

Noise  measurements  made  on  three  typical  NSVR  geometries  in- 
dicated that  the  signal-to-noise  ratio  of  a vortex  resistor  operating  in 
the  laminar  flow  regime  is  better  than  900. 


2 . LITERATURE  SURVEY 

The  concept  of  the  NSVR  (or  vortex  diode)  was  originated  by 
Dieter  Thoma  , and  the  SVR  configuration  was  first  used  by  the  Lee 
Company2  in  the  design  of  their  "Lee  Visco  Jet."  No  analytical  studies 
(in  the  laminar  regime)  of  the  SVR  uave  been  reported  in  the  open  lit- 
erature. With  the  exception  of  reference  , wherein  Paul  has  obtained 
order -of -magnitude  estimates  for  the  NSVR  performance,  no  analytical 
studies  have  been  reported.  The  vortex  and  diffuser  chambers  (see 
Figure  4)  have  been  studied  separately  as  described  in  the  sections 
below. 


2.1  Vortex  Flow 

Vortex  chamber  studies  have  been  motivated  by  applications  in 
fluidics  and  in  meteorology.  Four  different  techniques  have  been  ur.ed 
in  the  past  with  varying  degrees  of  success  for  characterizing  and  pre- 
dicting the  vortex  flow  field.  These  techniques  are: 

1.  Similarity  solutions, 

2.  Momentum  integral  solutions, 

3.  Finite  difference  solutions, 

4.  Empirical  approaches. 


^Thoma,  D.  "Fluid  Lines,"  U.S.  Patent  No.  1,839,618,  Patented 
5 January  1932. 

^Technical  Hydraulic  Handbook,  The  Lee  Comp-ny,  Westbrook, 
Connecticut,  1971. 

3Paul,  F.  W.  "Steady-State  Performance  of  Flueric  Diodes." 
(Ph.D.  dissertation,  Lehigh  University,  1968.) 
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Similarity  solutions  for  a potential  vortex  flow  over  a station- 
ary flat  plate  have  been  obtained  by  Kidd  and  Farris  , and  Hoffmann5. 

A numerical  search  is  required  to  solve  the  two-point  boundary  value 
problem  resulting  from  the  similarity  solution  approach.  It  has  been 
demonstrated  in  references  4 and  5 that  the  convergence  time  becomes 
too  large  for  tangential  Reynolds  numbers  greater  than  4.75.  Tan- 
gential Reynolds  numbers  for  the  SVR  and  the  NSVR  are  expected  to  be  of 
the  order  of  a few  hundreds  at  normal  operating  conditions.  Hence  a 
similarity  solution  technique  in  the  existing  form  does  not  seem  feasible. 


Von  Kantian' s integral  technique  (commonly  called  the  momentum 
integral  technique)  is  an  approximate  approach.  In  this  approach,  the 
partial  differential  equations  describing  the  flow  field  are  reduced  to 
ordinary  differential  equations  by  integrating  along  one  of  the  two 
independent  space  variables.  This  integration  (commonly  performed  along 
the  axis  of  the  vortex)  requires  velocity  profiles  in  functional  forms. 

Two  different  approaches  have  been  employed  to  obtain  reasonable  velocity 
profile  functions.  The  first  approach  involves  the  concept  of  Prandtl's 
boundary  layer.  Velocity  profile  functions  are  developed  or.  the  basis 
of  physical  insight  into  the  problem,  as  was  originally  done  by  Von  Karman. 
The  second  approach  due  to  Boyak  and  Rice^,  assumes  a general  polynominal 
form  for  the  velocity  profiles  without  requiring  any  a priori  knowledge 
of  the  flow  field.  In  this  respect,  the  second  approach  is  superior  to 
the  first  one  but  it  tends  to  be  computationally  inferior  due  to  the 
larger  number  of  differential  equations  involved  in  the  model. 


Kidd,  G.  J.,  Jr.  and  G.  J.  Farris.  "Potential  Vortex  Flow 
adjacent  to  a Stationary  Surface,"  ASME  Paper  No.  68-APM-15. 

^Hoffmann,  J.  A.  "Similarity  Solutions  for  the  Interaction  of 
a Potential  Vortex  with  Free  Stream  Sink  Flow  and  a Stationary  Surface," 
Trans . ASME,  J.  Fluids  Engineering,  Vol.  9i,  No.  1 (March,  1974), 
pp.  49-54. 

°Von  Karman,  T.  "Uber  laminar e and  turbulente  Reibung,"  Zeits , 
Ang.  Math.  Mech.  I (August,  1921),  pp.  233-252. 

^Boyack,  B.  E.,  and  W.  Rice.  "Integral  Method  for  Flow  between 
Rotating  Disks,"  Trans.  ASME,  J.  Basic  Engineering,  Series  D,  Vol.  93,  No.  3 
(September,  1971),  pp.  350-354^ 
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The  analyses  of  Rott  and  Lewellen®,  Wormley^,  and  Kwok,  Thinh 
and  Lin^O,  are  based  on  the  first  momentum  integral  approach.  Rott  and 
Lewellan's  model  (reference  8)  is  a general  one  applicable  to  both 
laminar  and  turbulent  flow  regimes.  Zero  radial  velocity  is  assumed  in 
the  core  (i.e.,  the  space  between  the  top  and  bottom  boundary  layers 
of  the  vortex  chamber)  in  order  to  produce  a simplified  model  amenable 
to  a closed-form  analytical  solution.  Wormleyy  removed  this  restriction 
by  including  a developing  flow  region  in  which  the  radial  velocity  in  the 
core  is  allowed  to  be  different  from  zero.  Wormley's  model  is  valid  for 
operation  in  the  turbulent  regime;  it  cannot  be  used  directly  for  model- 
ing the  SVR  and  the  NSVR  for  operation  in  the  laminar  regime.  But 
Wormley's  basic  approach  can  be  utilized  by  including  laminar  shear 
stresses  and  velocity  profile  functions  which  are  appropriate  for  the 
laminar  regime. 

A more  approximate  form  of  the  momentum  integral  approach,  as 
applied  to  vortex  flow,  is  the  "one-dimensional"  approach.  In  this 
approach,  the  velocity  profiles  are  assumed  to  be  uniform.  Closed  form 
solutions  based  on  the  one-dimensional  approach  have  been  obtained  by 
Taplin^,  and  Bichara  and  Ornerl^.  Taplin^  treated  the  tangential 
velocity  as  a known  fuction  of  the  radial  distance.  Bichara  and  Orner^ 
developed  a more  general  model  (turbulent  regime)  by  calculating  the 
radial  distribution  of  the  velocities  using  appropriate  equations.  This 
latter  approach  is  highly  efficient,  and  can  be  extended  for  use  in  the 
laminar  regime.  But,  the  results  should  be  expected  to  have  limited 
validity  due  to  the  uniform  velocity  profile  assumption. 

The  availability  of  high  speed  electronic  computers  with  rela- 
tively large  storage  capacities  makes  finite-difference  methods  of  solv- 
ing partial  differential  equations  more  attractive.  A recent  application 
of  this  approach  related  to  the  vortex  flow  problem  was  presented  by 


8Rott,  N.,  and  W.  S.  Lewellen.  "Boundary  Layers  and  Their  Inter- 
actions in  Rotating  Flows,"  Progress  in  Aeronautical  Sciences,  Vol.  7 
(1966),  pp.  111-144. 

^Wormley,  D.  N.  "Analytical  Model  for  the  Incompressible  Flow 
in  Short  Vortex  Chambers,"  Trans . ASME,  J.  Basic  Engineering,  Series  D, 

Vol.  91,  No.  1 (June,  1969),  pp.  2bU-2lC. 

10Kwok,  C.  C.  K.,  N.  D.  Thinh,  and  S.  Lin.  "An  Investigation  of 
Confined  Vortex  Flow  Phenomena,"  Trans . ASME,  J,  Basic  Engineering, 

Series  D,  Vol.  94,  No.  3 (September,  1972),  pp.  689-696. 

^Taplin,  L.  B.  "Small  Signal  Analysis  of  Vortex  Amplifiers," 
Chapter  7,  AGARDograph,  118  (December,  1968),  pp.  235-295. 

12gichara,  R.  T. , and  P.  A.  Orner.  "Analysis  and  Modeling 
of  the  Vortex  Amplifier,"  Trans . ASME, J.  Basic  Engineering,  Series  D, 

Vol.  91,  No.  1 (December,  1969),  pp.  755-763. 
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Farris  et  al.  The  computational  times  reported  in  reference  indi- 
cate that  the  state-of-the-art  in  this  area  of  numerical  computation  has 
not  yet  reached  a reasonable  level  of  practicality. 

Very  few  empirical  or  "black-box"  studies  of  the  vortex  flow 
problem  have  been  reported  in  the  literature.  The  most  relevant  here  is 
the  empirical  study  of  the  vortex  amplifier  by  Hondour  and  Rezek^.  The 
important  effect  of  viscosity  variation  of  the  hydraulic  oil  (due  to 
temperature  changes)  was  not  included  in  this  model  however.  Since  the 
vortex  amplifier  characteristics  are  strongly  influenced  by  viscosity 
variations,  the  model  presented  in  reference  14  is  of  limited  use. 

Studies  conducted  on  vortex  rate  sensors  are  not  pertinent  to 
the  study  of  the  SVR  or  the  NSVR  due  to  the  significant  differences  in 
the  associated  flow  fields  observed  in  these  two  classes  of  devices.  The 
Rossby  number  (ratio  of  radial  to  tangential  velocity)  is  very  large 
(much  greater  than  unity)  for  a vortex  rate  sensor  and  it  is  of  the  order 
of  0.1  for  the  SVR  and  the  NSVR. 

In  summary,  no  efficient  analysis  technique  has  been  reported 
which  is  directly  applicable  to  the  vortex  chamber  under  normal  opera- 
ting conditions.  An  empirical  study  of  the  vortex  chamber  is  feasible. 

An  alternative  semi-analytical  approach  is  the  one-dimensional  approach 
of  Bichara  and  Orner  with  necessary  modifications  for  operation  in 
laminar  regime.  Yet  another,  but  more  difficult,  approach  is  to  adapt 
Wormley's  momentum  integral  approach  to  the  laminar  flow  case. 


2.2  Orifice  Flow 

The  axial  flow  through  the  central  orifice  in  the  vortex  resis- 
tor is  accompanied  by  a large  swirl  at  the  inlet  induced  by  the  flow  in 
the  vortex  chamber  (forward  flow  case).  This  swirl  flow  through  the 
orifice  resembles  the  swirl  flow  through  nozzles  used  in  swirl  atomizers. 
The  effect  of  swirl  is  to  force  a large  portion  of  the  fluid  to  the 
boundary  layer  on  the  cylindrical  surface,  and  thereby  increase  the 
resistance  to  flow  through  the  orifice.  In  the  NSVR  there  is  no  swirl 
in  the  orifice  when  flow  is  in  the  "reverse"  direction  (i.e.,  low  resis- 
tance direction). 


13Farris , G.  J.,  G.  J.  Kidd,  Jr.,  D.  W.  Lick,  and  R.  E.  Textor. 
"A  Theoretical  and  Experimental  Study  of  Confined  Vortex  Flow,"  Trans . 
ASME,  J.  Appl.  Mech.,  Series  E,  Vol,  91,  No.  4 (December,  1969), 
pp.  687-692. 

14 Hondour , Y.  T. , and  C.  W.  Rezek.  "Hydraulic  Vortex  Valve, 
Modeling  and  Design,"  ASME  Paper  No.  74-WA/Aut- 19 , (1974). 
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Taylor^  and  Weber^  have  studied  the  swirl  flow  through  noz- 
zles using  the  momentum  integral  approach.  It  is  possible  to  adapt  this 
approach  for  analyzing  the  central  orifice  in  the  vortex  resistor.  But, 
the  resulting  model  would  significantly  increase  the  complexity  of  the 
overall  model  of  the  vortex  resistor.  A simple  model  (essentially  one 
dimensional)  which  takes  into  account  the  effect  of  swirl  in  the  orifice 
has  been  developed  by  Bichara  and  Orner^.  Since  the  inlet  conditions 
for  the  orifice  are  not  exactly  the  ideal  conditions  used  in  the  more 
sophisticated  approach  of  references  15  and  16,  it  appears  more  reason- 
able to  use  the  approximate  but  efficient  model  of  reference  12. 

The  vortex  resistor  rectangular  cross-section  orifices  (i.e., 
the  inlet  orifice  for  the  NSVR  and  the  inlet  and  outlet  orifices  for 
the  SVR)  and  circular  cross-section  central  orifice,  can  be  modeled  as 
short-tube  orifices  with  sharp-edged  entries.  The  coefficient  of  dis- 
charge reported  by  Lichtarowicz,  et  al.^  and  Reid  can  be  used  for 
predicting  the  flow  characteristics  of  these  orifices. 


2.3  Diffuser  Flow 


The  diffuser  chamber  in  the  case  of  the  NSVR  is  a sudden  expan- 
sion section.  This  chamber  can  be  modeled  with  reasonable  accuracy  by 
means  of  a one-dimensional  momentum  equation.  The  diffuser  chamber  of  the 
SVR  is  characterized  as  a decelerating  vortex  flow  (counter-rotating  with 
respect  to  the  inlet  vortex  chamber);  the  corresponding  flow  field  is 
exceedingly  complex. 


12Bichara,  R.  T. , and  P.  A.  Orner.  "Analysis  and  Modeling  of 
the  Vortex  Amplifier,"  Trans . ASME,  J- Basic  Engineering,  Series  D, 

Vol . 91,  No.  1 (December,  1969),  pp.  7c"'-"63. 

15Taylor,  G.  I.  "The  Boundary  Layer  in  the  Converging  Nozzle  of 
a Swirl  Atomizer,"  Quarterly  J.  Mechanics  and  Applied  Mathematics, 

Vol.  3 (1950),  pp.  129-140. 

16Weber,  H.  E.  "The  Boundary  Layer  Inside  a Conical  Surface  due 
to  Swirl,"  Trans . ASME,  J.  Applied  Mechanics,  Series  E,  Vol.  78,  No.  4 
(December,  1956),  pp.  587-592. 

^ 'Lichtarowicz , A.,  R.  K.  Duggins,  and  E.  Markland.  "Discharge 
Coefficients  for  In-compressible  Non-Cavitating  Flow  Through  Long  Ori- 
fices," J.  Mechanical  Engineering  Sciences,  Vol.  7,  No.  2 (1965),  pp.  210- 
219. 

^Reid,  K.  N.  "Static  Characteristics  of  Fluid  Amplifiers," 
Proceedings  of  the  Fluid  Power  Research  Conference,  Oklahoma  State 
University,  July  25-26,  1967. 
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Radial  diffuser  flows  have  been  studied  both  analytically  and 
experimentally  by  Jansen  19,20.  His  analyses  are  based  on  a momentum 
integral  technique.  The  flow  field  in  the  radial  diffuser  is  assumed 
to  be  symmetric  about  a mid-plane,  with  good  accuracy.  But,  the  flow 
field  in  the  diffuser  chamber  of  the  SVR  is  dominantly  anti-symmetric 
and  resembles  the  flow  field  on  one  side  of  an  enclosed  rotating  disc. 
Therefore,  Jansen's  results  cannot  be  used  for  the  vortex  resistor 
case . 


There  is  no  published  literature  on  the  analysis  of  the  flow 
field  in  a diffuser  chamber  such  as  that  in  an  SVR.  But,  the  rotating 
disc  studies  can  serve  as  a guide.  Analytical  and  experimental  studies 
of  an  enclosed  disc  were  originally  conducted  by  Daily  and  Nece2!. 
Recently  Miloh  and  Poreh22  have  extended  the  analysis  to  non-Newtonian 
fluids.  The  analysis  involves  consideration  of  two  boundary  layers,  one 
on  the  rotating  disc  and  the  other  on  the  stationary  casing.  Rigid  body 
rotation  is  assumed  for  the  core  and  a momentum  integral  approach  is 
used  to  solve  for  the  flow  field.  These  analyses  do  not  include  through- 
flow.  The  basic  approach  of  Daily  and  Nece  can  be  extended  to  the 
diffuser  chamber  of  the  vortex  resistor.  Alternatively,  one -dimensional 
approach,  similar  to  the  one  of  reference  12,  can  be  used  to  approximate 
the  complex  flow  field  existing  in  the  diffuser  chamber  of  the  SVR. 


2 .4  Vortex  Dynamics 

No  analytical  or  experimental  results  have  been  reported  in  the 
open  literature  concerning  the  dynamical  behavior  of  SVR  and  NSVR  elements. 


^Bichara,  r.  T. , and  P.  A.  Orner.  "Analysis  and  Modeling  of 
the  Vortex  Amplifier,"  Trans.  ASME , J.  Basic  Engineering,  Series  D,  Vol. 

91,  No.  1 {December,  1969),  pp.  755-763. 

™Jansen,  w.  "Rotating  Stall  in  a Radial  Vaneless  Diffuser," 
Trans . ASME , J.  of  Basic  Engineering,  Series  P,  Vol.  86,  No.  4 (December, 
1964),  pp.  750-758. 

uJansen,  W.  "Steady  Fluid  Flow  in  a Radial  Vaneless  Diffuser," 
Trans . ASME , J.  Basic  Engineering.,  Series  D,  Vol.  86,  No.  3 (September, 
1964) , pp.  607-619. 

2lDaily,  J.  W. , and  R.  E.  Nece.  "Chamber  Dimension  Effects  on 
Induced  Flow  and  Frictional  Resistance  of  Enclosed  Rotating  Disks,”  Trans . 
ASME,  J Basic  Engineering,  Series  D,  Vol.  82,  No.  1 (March,  1960), 
pp"I  217-232 . 

22Miloh,  T.,  and  M.  Poreh.  "The  Resistance  to  Rotation  of  Free 
Enclosed  Disks,"  Trans . ASME,  J. Appl.  Mech. , Series  E , Vol.  93,  No  4 
(December,  1971),  pp.  749-755. 
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However,  some  work  has  been  done  on  the  dynamics  of  fluidic  vortex 
rmplifier  and  vortex  rate  sensors. 

The  vortex  rate  sensor  has  a low  swirl  component  whereas  a 
NSVR  inherently  has  a large  swirl  component.  Results  of  dynamic  studies 
of  vortex  rate  sensors  23,  24,  and  25  appear  to  be  of  limited  value  in 
understanding  the  dynamical  behavior  of  vortex  resistance  elements. 
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Duff,  et  al.  and  Knapp  have  measured  the  dynamic  response 
of  a vortex  valve,  but  no  attempt  has  been  made  to  analytically  predict 
the  response.  Taplin^  has  derived  an  analytical  model  for  small-signal 
operation  of  a vortex  amplifier.  His  "one-dimensional  model"  includes 
a time  delay  Tj  and  a first  order  lag  ?2  to  characterize  the  vortex 
chamber  dynamics.  The  propagation  of  tangential  velocity  from  the  outer 
edge  of  the  vortex  chamber  to  the  exit  orifice  gives  rise  to  the  pure 
time  delay  T^.  For  inviscid  flow,  T^  equals  the  chamber  fill  time.  But, 
due  to  viscous  effects,  the  flow  is  contained  largely  in  the  boundary  layers 
on  the  top  and  bottom  surfaces  of  vortex  chamber.  This  tends  to  reduce 
the  effective  chamber  volume  (see  reference  29).  Taplinz  used  one 
quarter  of  the  chamber  fill-time  as  the  time  delay.  The  second  time 
constant  T2  arises  from  the  inertia  of  the  core  between  the  two  boundary 
layers.  Again  one  quarter  of  the  chamber  fill-time  was  used  by  Taplin^o 
for  T2. 


^Sarpkaya,  T.  "A  Theoretical  and  Experimental  Investigation  of 
the  Vortex-Sink  Angular  Rate  Sensor,"  Proceedings  of  the  HDD  Amplification 
Symposium,  Vol.  II  (October,  1965). 

2^Sarpkaya,  T. , J.  M.  Goto,  and  J.  M.  Kirshner.  "A  Theoretical 
and  Experimental  Study  of  Vortex  Rate  Gyro,"  Advances  in  Fluidics,  (May,  1967). 

^ostdiek,  A.  J.  Viscous  Vortex  Rate  Sensor,  Technical  Report  No. 
TR-1555,  Harry  Diamond  Laboratories,  (November,  1971). 

2^Duff,  J.,  K.  Foster,  and  D.  G.  Mitchell.  "Some  Experiments 
on  the  Vortex  Valve,"  First  Cranfield  Fluidics  Conference,  British  Hydro- 
mechanical Research  Association,  Cranfield,  England  (1965). 

2?Knapp,  A.  B.  "The  Development  and  Control  of  a Vented  Vortex 
Amplifier,"  Third  Cranfield  Fluidics  Conference,  British  Hydromechahical 
Research  Association,  Cranfield,  England,  (1968). 

^^Taplin,  L.  B.  "Phenomenology  of  Vortex  Flow  and  its  Appli- 
cations to  Signal  Amplification,"  Paper  at  Engineering  Seminar  on  Fluid 
Control  Systems,  Pennsylvania  State  University  Engineering  Proceedings 
(1967),  p.  45. 

29xirshner,  J.  M. , and  S.  Katz.  Design  Theory  of  Fluidic 
Components . New  York:  Academic  Press,  1975,  pp.  276-310. 
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In  contrast,  Bell®^  used  T^  = 1/2  the  chamber  fill-time,  and  T2  = 1.3 
times  the  chamber  fill-time  to  match  his  experimental  data.  The  empirical 
nature  of  the  models  used  by  Taplin  and  Bell  is  obvious. 

To  rationalize  the  choices  of  T^  and  T£,  Anderson^  derived  a 
distributed  parameter  dynamic  model.  The  model  was  used  to  estimate  the 
times  T^  and  T2 . The  resulting  times  were  fractions  of  the  chamber  fill- 
time ranging  from  0.2778  to  0.3648  for  T^  and  from  0.6454  to  0.7647  for 
T2  for  the  four  cases  studied  by  Anderson.  Unlike  Bell^O  and  Taplin28j 
Anderson  does  not  predict  a constant  fraction  of  the  chamber  fill-time 
for  Ti  and  T2.  In  this  respect,  Anderson's  model  is  superior  to  the  others. 
Experimental  results  confirm  the  validity  of  Anderson's  approach. 

Anderson's  model  is  based  in  part  on  Wormley's  steady-state 
modely  which  in  turn  uses  turbulent  shear  stress  terms.  A major  modifica- 
tion would  be  required  in  the  Anderson  model  to  make  it  useful  in  the 
dynamic  characterization  of  the  NSVR  and  SVR  which  typically  operate  in 
the  laminar  flow  regime. 

Although  the  previous  dynamic  studies  have  limited  direct  appli- 
cation to  the  vortex  resistor  problem,  Taplin^®  and  Bell30  suggest  a model 
form  for  describing  the  vortex  chamber  dynamics.  That  is,  a model  form 
which  includes  a pure  time  delay  and  a first  order  lag  agrees  reasonably 
well  with  experimental  data. 


3.  STEADY-STATE  EMPIRICAL  MODELS 


In  cases  where  fluid  systems  involving  a large  number  of  elements 
are  to  be  simulated,  it  is  often  desirable  to  use  simple  mathematical  models 
for  the  elements  in  order  to  maximize  computational  efficiency.  The  empir- 
ical models  developed  here  for  symmetric  and  non- symmetric  vortex  resistors 
are  single  algebraic  equations  which  contain  experimentally  derived  coefficients. 
They  are  computationally  efficient  for  preliminary  design  calculations  but 
limited  in  their  range  of  application  by  the  available  data  base.  More 
accurate  analytical  models  are  developed  in  section  4 which  consist  of  sets 


^Wormley,  D.  N.  "Analytical  Model  for  the  Incompressible  Flow  in 
Short  Vortex  Chambers,"  Trans.  ASME,  J.  Basic  Engineering,  Series  D,  Vol.  91, 
No.  1 (June,  1969),  pp.  264-276. 

28Taplin,  L.  B.  "Phenomenology  of  Vortex  Flow  and  its  Applications 
to  Signal  Amplification,"  Paper  at  Engineering  Seminar  on  Fluid  Control 
Systems,  Pennsylvania  State  University  Engineering  Proceedings  (1967),  p.  45. 

®°Bell,  A.  C.  "Optimization  of  a Vortex  Valve," (M.S.  thesis, 
Massachusetts  Institute  of  Technology,  1965.) 

-^Anderson,  W.  W.  "A  Dynamic  Model  of  Vortex-Type  Fluid  Amplifier," 
(Ph.D.  dissertation,  Massachusetts  Institute  of  Technology,  1971.) 
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of  algebraic/differential  equations;  these  models  are  applicable  over  a wide 
range  of  geometrical  and  operational  conditions,  but  they  require  the  use  of 
a digital  computer  in  performing  analysis  and  design  calculations. 


3.1  Model  Form 

The  empirical  models  developed  in  this  study  are  similar  in  func- 
tional form  to  the  vortex  amplifier  model  developed  by  Hondour  and  Rezekl^ 
except  that  the  effect  of  viscosity  is  included  here.  The  steady-state  char- 
acteristics of  a vortex  resistor  is  assumed  to  be  of  the  form 
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where  the  a's  are  constants  which  are  computed  by  linear  regression  using 
experimental  data.  The  form  of  equation  (1)  is  consistent  with  that  normally 
used  in  analyses  and  data  correlations  for  orifice  type  resistors. 


3.2.  Experimental  Apparatus  and  Procedure  for  Measurement 


of  Steady-State  Characteristics 


Figure  9 shows  a schematic  of  the  experimental  apparatus  used  to 
obtain  steady-state  characteristics  for  vortex  resistors  over  a range  of 
geometrical  and  operational  conditions.  These  data  were  used  to  determine 
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Figure  9.  Schematic  of  experimental  apparatus  for  measurement  of 
steady-state  characteristics. 


i4Hondour,Y.  T. , and  C.  W.  Rezek.  "Hydraulic  Vortex  Valve, 
Modeling  and  Design,"  ASME  Paper  No.  74-WA/Aut-19 , (19?4). 
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the  constants  in  equation  (1)  and  to  validate  the  one -dimensional  and  two- 
dimensional  models  presented  in  section  4.  A Vickers  hydraulic  piston 
pump  rated  at  6.3  m^/s  (at  1.03x10^  kN/m2)  provided  the  fluid  power  source. 
The  operating  fluid  was  MIL-5606A  hydraulic  oil.  A chromalox  coil- type 
immersion  hefter  (of  10  kW  capacity)  was  used  for  heating  the  hydraulic 
oil  (35* -57.2  C).  A feedback  temperature  controller  was  used  to  maintain 
the  oil  temperature  constant  at  any  desired  value  above  the  room  temperature 
Figure  10  is  a photograph  of  the  test  stand  and  primary  instrumentation. 


Figure  10.  Photograph  of  the  test  stand  and  primary  instrumentation. 


Resistor  Block 


The  vortex  resistor  block  was  of  laminated  construction.  Figure 
11  shows  parts  of  the  NSVR  assembly.  Tapped  holes  in  the  cop  and  bottom 
blocks  allowed  mounting  of  the  pressure  gauges  and  thermocouples  very  near 
the  effective  inlet  and  exit  of  the  resistor.  Locating-pins  were  used  in 
accurately  drilled  holes  to  insure  proper  alignment  of  the  laminates. 

Figure  11  shows  the  vortex  chamber  orifice  and  diffuser  laminations  used 
for  the  SVR  configuration.  (The  vortex  chamber  and  diffuser  chamber  lam- 
inations are  identical.)  For  the  non-cymmetric  configuration  or  NSVR,  no 
diffuser  laminations  were  used.  The  end  block  of  the  NSVR  contained  a 
sudden  expansion  diffuser  aligned  with  the  central  orifice  as  shown  in  Figure 
12.  Figure  13  is  a photograph  showing  the  various  parts  of  the  NSVR  assembly. 
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Laminated  construction  for  the  SVR 
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Figure  13.  Exploded  view  of  the  NSVR  block  assembly 


Table  I lists  the  various  sets  of  laminations  used  in  this  study 
Sets  1 and  2 were  made  by  Air-Research  Corporation.  Sets  3 and  4 were 
designed  and  constructed  at  OSU.  Sets  1 and  2 were  made  from  steel  and 
were  copper  plated.  Set  3 was  made  of  brass  and  set  4 of  plexi-glass. 

In  sets  1 through  3,  surfaces  in  contact  were  sufficiently  smooth  and 
the  clamping  force  during  assembly  was  large  enough  to  result  in  negli- 
gible leakage  between  the  laminations.  Rubber  0-ring  seals  were  used 
with  set  4 laminations  to  eliminate  leakage.  Different  end  blocks  were 
used  for  the  four  sets  of  laminations  due  to  their  different  geometric 
scales.  Different  overall  resistor  geometries  were  assembled  from  sets 
1,  2,  and  3.  The  height  of  the  vortex  chamber  (and  diffuser  chamber  in 
the  case  of  the  SVR)  and  the  length  of  the  central  orifice  were  conven- 
iently varied  by  the  number  of  laminations  used. 


Instrumentation 


The  following  quantities  were  measured 


Inlet  and  outlet  static  pressures, 
Volumetric  flow  rate  through  the  resistor 
Temperature  of  the  hydraulic  oil. 


Bourdon- tube  type  pressure  gauges  were  used  to  measure  the  static 
pressures  in  the  range  3.5  kN/m^  and  up.  For  pressures  below  3.5  kN/m  , 
strain-gauge  pressure  transducers  were  used  for  improved  accuracy.  The 
pressure  gauges  were  calibrated  with  a dead-weight  tester.  The  maximum 
error  in  pressure  measurement  was  + 0.2%  of  the  full  scale. 

Fischer-Porter  rotameter-type  flowmeters  were  used  to  measure 
the  volumetric  flow  rate.  Three  different  flowmeters  were  used  to  cover 
the  entire  range  of  flow  rate,  i.e. , 6x10"®  to  6.5xl0"5  m^/s.  The  flow- 
meters were  calibrated  with  the  actual  operating  fluid  over  the  entire 
range  of  temperature  operation.  Maximum  error  in  flow  rate  measurement  was 
0.47,  of  the  full  scale. 

A Leeds  and  Northrup  Numatron  II  temperature  measurement  system 
was  used  in  conjunction  with  locally  built  Iron-Constantan  thermocouples. 
The  thermocouple  elements  were  located  in  the  end  blocks  of  the  vortex 
resistor.  The  flow  sensitivities  of  the  thermocouples  were  minimal  (of 
the  order  of  0.5*C  error  for  the  total  range  of  flow  rate  used)  when  they 
were  located  in  the  end  block  position  shown  in  Figure  14.  The  overall 
error  of  the  temperature  measurement  system  was  less  than  t 0.5*C. 


Figure  14.  Location  of  the  thermocouple  in  the  end  block. 

The  temperature  measurement  system  was  calibrated  by  measuring 
the  temperature  of  boiling  water.  Two  thermocouples,  one  at  the  upstream 
and  one  at  the  downstream  of  the  resistor,  were  used  to  determine  the 
temperature  difference  across  the  resistor.  By  thermally  insulating  the 
resistor  block  the  temperature  difference  in  the  worst  case  (at  lowest 
flow  rate)  was  t 1.4*c.  An  average  of  the  two  temperatures  was  used  as 
the  temperature  of  the  hydraulic  oil  for  computational  purposes. 
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Density  and  viscosity  were  the  only  properties  of  the  hydraulic 
oil  used  in  the  analysis  of  the  resistor.  Properties  data  available  from 
the  literature  (see  Figures  15  and  16)  were  used  for  all  calculations. 

The  pressure  dependence  of  viscosity  and  density  was  not  included  due  to 
the  relatively  low  operating  pressure.  Care  was  taken  to  remove  all 
trapped  air  from  the  system  by  bleeding. 


Specific 

Gravity 


Figure  15.  Specific  gravity  of  MIL-5606A  as  a function  of  temperature. 

(See  reference  2) 
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Figure  16.  Kinematic  viscosity  of  MIL-5606A  as  a function  of 
temperature.  (See  reference  2) 


All  measured  steady-state  data  are  presented  in  detail  in  a 
supplementary  report*  to  HDL.  Typical  pressure  drop-flow  rate  character- 
istics are  shown  in  section  4. 


3.3.  Empirical  Correlations 

Appendix  A outlines  the  linear  regression  procedure  used  to 
obtain  the  model  constants  (a^'s  of  Equation  1 based  on  the  measured 
data).  Figure  17  shows  a correlation  for  a NSVR  with  a radius  ratio 
(^/r^)  of  0.299.  Other  parameters  for  this  case  are  listed  in  Table  III 
in  Appendix  A.  The  flow  conditions  ranged  from  = 100  to  1000.  The 

model  agrees  with  the  experimental  data  within  30  percent  of  reading. 

X102 


Figure  17.  Empirical  model  of  a NSVR  for  operation  in  the 
high  resistance  flow  direction. 


4.  STEADY-STATE  ANALYTICAL  MODELS 

This  chapter  summarizes  the  development  of  approximate  analytical 
models  which  predict  the  steady-state  behavior  of  the  two  types  of  vortex 
resistors.  The  approach  is  to  approximate  the  basic  flow  processes  with 


*Steady-State  Characteristics  of  Vortex  Resistors:  "Experimental 
Data,"  Supplement  to  the  report  submitted  by  Oklahoma  State  University  to 
Harry  Diamond  Laboratories  under  contract  No.  DAAG-39-75-C-0028 , February,  1976. 


simple  and  manageable  functions.  In  the  case  of  the  vortex  chamber  (see 
Figure  4) , the  flow  field  has  been  relatively  well  described  in  earlier 
studies^.  The  flow  field  in  the  diffuser  chamber  of  the  SVR  has  not 
been  studied  and  is  difficult  to  visualize  intuitively.  A flow  visuali- 
zation study  was  conducted  to  identify  the  basic  character  of  the  flow 
field  in  the  symmetric  vortex  resistor. 


4.1.  Flow  Visualization  Stud 


A large  scale  plexi-glass  model  (50.8mm  outer  diameter  vortex 
chamber)  run  on  water  was  used  for  the  flow  visualization  study.  Figure 
18  shows  the  schematic  of  the  experimental  setup.  Figure  19  shows  an 
exploded  view  of  the  model  with  the  dye-injection  system.  The  first 
plate  on  the  left  carries  an  injector  (inside  the  tube)  which  was  used 
for  introducing  dye  at  different  depths  in  the  vortex  chamber.  The 
vortex  chamber,  central  orifice  and  diffuser  plates  can  also  be  seen  in 
Figure  19.  The  flow  conditions  were  maintained  at  a Reynolds  number 
(based  on  the  central  orifice  diameter)  of  200  for  all  cases  except  that 
shown  in  Figure  27.  For  Figure  27,  the  Reynolds  number  was  500. 
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Figure  18.  Schematic  of  flow  visualization  apparatus. 


^Wormley,  D.  N.  "Analytical  Model  for  the  Incompressible  Flow 
in  Short  Vortex  Chambers,"  Trans.  ASME,  J.  Basic  Engineering,  Series  D. 
Vol.  91,  No.  1 (June,  1969),  pp.  264-276. 


Figure  19.  Plexiglass  model  for  flow  visualization 


4.1.1.  Vortex  Chamber 

Figures  20-23  show  typical  results  from  the  flow  visualization 
study  of  the  vortex  chamber.  Figures  20  and  21  show  the  boundary  layer  flow. 


$ 


Figure  20.  Boundary  layer  flow  in  vortex  chamber. 

In  the  cases  shown  in  Figures  22  and  23,  dye  was  introduced  into  the  core 
region.  The  dye  streaks  in  these  latter  cases  show  the  core  flow  with  a 
relatively  large  tangential  to  radial  velocity  ratio.  The  streamlines 
appear  to  be  disturbed  due  to  the  wake  created  by  the  injector. 
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Core  flow  in  vortex  chamber 


Diffuser  Chamber 


Figures  24-29  show  the  flow  patterns  in  the  diffuser  chamber. 
In  all  cases  the  diffuser  chamber  was  first  filled  with  dye,  and  then 
water  flow  was  initiated  by  opening  the  needle  valve  (see  Figure  18). 

The  flow  patterns  of  Figures  24,  25,  and  26  were  observed  in  time  succes 
sion.  A spinning,  tube-like  flow  pattern  was  observed  during  the  test. 
(The  spinning  of  the  tube  was  about  a toroidal  axis  [see  Figure  30] ) . 


Figure  27.  Flow  pattern  in  the  diffuser  chamber  of  SVR 
at  a Reynolds  Number  of  500. 


Figure  28.  Reverse  flow  in  the  diffuser  chamber  of  SVR 
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Figure  29  shows  the  boundary  layer  flow  along  the  annular  wall  and  reverse 
flow  on  the  top  surface  of  the  diffuser  chamber.  An  attempt  to  photo- 
graph cases  where  dye  was  injected  into  the  core  region  of  the  diffuser 
failed.  The  dye  injected  at  various  heights  dispersed  quickly  without  a 
stable  pattern.  The  presence  of  the  injector  and  the  rate  of  injection 
disturbs  the  flow  field  strongly.  Recirculation  was  observed  in  the  core 
region  of  the  orifice.  Figure  30  shows  schematically  the  overall  flow  field 
observed  in  the  SVR. 


Figure  29.  Reverse  flow  in  the  diffuser  chamber  of  SVR. 


The  results  shown  in  Figures  20-23  indicate  that  a boundary  layer- 
vortex  core  model  would  be  an  appropriate  description  of  the  flow  field 
in  the  vortex  chamber.  Figure  22  illustrates  that  the  flow  in  the  diffuser 
chamber  is  not  symmetric  about  a mid-plane.  The  conceptual  model  for  the 
diffuser  chamber  ideally  should  account  for  the  observed  reverse  flow. 
Boundary  layer  flows  may  be  assumed  on  the  top  and  bottom  surfaces  of  the 
diffuser  chamber  with  the  radial  component  of  velocity  being  radially 
inward  at  the  top  and  radially  outward  at  the  bottom. 
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Figure  30.  Flow  patterns  in  vortex  resistor. 


One-Dimensional  Model  of  the  NSVR 


Model  Description 


The  goemetry  of  the  NSVR  may  be  divided  into  four  regions  follow- 
ing the  approach  of  Bichara  and  Orner*2.  Figure  31  shows  the  four  regions. 
These  regions  may  be  modeled  by  making  assumptions  appropriate  for  the 
individual  regions.  The  resulting  equations  may  be  suitably  combined  to 
result  in  an  overall  mathematical  model. 

The  modeling  approach  for  the  NSVR  is  l immarized  below  (see 
Appendix  B for  additional  d etails) . The  inlet  region  of  the  NSVR  is 
modeled  as  an  orifice.  The  discharge  coefficient  characteristic  is 


12Bichara,  R.  T.,  and  P.  A.  Orner.  "Analysis  and  Modeling  of  the 
Vortex  Amplifier,"  Trans . ASME.  J.  Basic  Engineering,  Series  D,  Vol.  91, 

No.  1 (December,  1969),  pp.  755-763. 
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Figure  31.  Regions  of  NSVR. 

assumed  to  be  the  same  as  that  of  a short-tube  orifice  with  a sharp-edged 
entry.  Experimental  data  on  the  circular  short-tube  orifice  reported  in 
reference  [18]  are  used  even  though  the  downstream  conditions  are  not  the 
same  in  the  two  cases.  An  empirical  correlation  (see  Appendix  B,  Equation 
(B.5))  was  developed  to  conveniently  use  the  experimental  orifice  data. 

The  flow  is  assumed  to  be  axisymmetric  at  the  periphery 
of  the  region  2.  The  shear  stress  on  the  outer  cylindrical  surface  of 
the  vortex  chamber  is  obtained  from  an  analogy  with  the  flat  plate  solu- 
tion of  Blassius-^.  This  assumption  is  understood  to  be  a rough  approxi- 
mation. A more  exact  approach  would  take  into  account  the  non-axisymmetric 
nature  of  the  flow  field  at  the  entry  of  the  vortex  chamber. 

Flow  within  region  2 of  the  NSVR  is  assumed  to  be  one- 
dimensional, i.e.,  the  static  pressures  and  velocities  are  assumed  to  be 
functions  of  the  radial  variable  "r"  only.  The  ordinary  differential 
equations  which  result  (see  Appendix  B)  are  much  simpler  than  the  corres- 
ponding ordinary  differential  equations  that  result  from  a momentum 


Schlichting,  H.  Boundary-Layer  Theory.  6th  ed.,  New 
York:  McGraw-Hill  Book  Co.,  Inc.,  1941. 
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integral  approach.  These  equations  can  be  solved  in  closed  form,  which 
eliminates  numerical  integration  and  results  in  high  computational  effi- 
ciency. The  assumption  of  uniform  velocity  profile  does  not  permit 
inclusion  of  the  shear  stress  into  the  model.  However,  a shear  stress 
term  can  be  artificially  included.  Radial  shear  is  neglected  but  shear 
stress  in  the  tangential  direction  is  included. 

The  presence  of  a developing  flow  region  in  the  vortex  chamber 
suggests  that  the  average  shear  stress  should  be  greater  than  the  shear 
stress  for  the  fully-developed  flow  case.  The  tangential  shear  stress 
value  was  treated  as  a "free  parameter,"  and  was  chosen  such  that  the 
computed  results  agreed  reasonably  well  with  experimental  data.  A value 
two  times  that  associated  with  fully-developed  laminar  flow  between  two 
parallel  plates  was  found  to  produce  acceptable  results. 

Region  3 is  modeled  as  an  orifice.  Separation  at  the  entry  is 
minimized  by  the  swirl,  and  the  wall  boundary  layer  is  thinned  due  to 
centrifugal  force  produced  by  the  swirl.  These  effects  tend  to  result 
in  a larger  orifice  discharge  coefficient  than  for  the  case  with  no  swirl. 
(For  the  case  with  no  swirl,  the  discharge  coefficient  varies  from  0.6 
to  0.8  for  the  range  of  operation  considered.)  A value  of  unity  was  used 
for  the  discharge  coefficient  of  the  central  orifice.  The  static  pressure 
at  the  downstream  end  of  the  central  orifice  is  assumed  to  be  constant 
over  the  area  of  the  orifice.  But,  the  upstream  pressure  is  allowed  to 
vary  with  the  radial  distance  r. 

Region  4 is  modeled  as  a sudden- expans ion  diffuser  which  recovers 
half  of  the  kinetic  energy  of  the  fluid.  This  choice  was  based  on  an 
approximate  analysis  of  the  diffuser  taking  into  account  the  effect  of 
swirl.  (See  Appendix  B for  details.) 


4.2.2.  Typical  Results 

Figures  32-35  compare  experimental  data  and  analytical  predictions 
based  on  the  one -dimensional  model.  The  model  predicts  the  experimentally 
observed  behavior  of  the  NSVR  reasonably  well  in  both  the  positive  and  neg- 
ative temperature  sensitivity  regions  for  a majority  of  the  geometries 
tested.  Theo temperature  range  for  the  analytical  predictions  was  limited 
to  35  to  57  C because  of  test  stand  limitations. 

Figures  36  & 37^  show  analytical  predictions  for  the  temperature 
range  -45.55  to  +65.55  C (a  common  military  specification).  The  geometry 
of  Figure  36  exhibits  a relatively  large  negative  temperature  sensitivity 
(-1.533x10*7  m-Vs-'c).  In  contrast,  Figure  37  shows  nearly  zero  tempera- 
ture sensitivity  (-5.900x10"^  m^/s-  C)  over  the  entire  temperature  range. 
Experimental  verification  is  needed  to  evaluate  the  accuracy  of  the  pre- 
dictions over  the  wide  temperature  range  used  in  Figures  36  and  37. 
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Figure  38  shows  the  steady-state  characteristic  of  the  inlet 
orifice  alone.  The  inlet  orifice  alone  always  exhibits  a positive  tem- 
perature sensitivity  & Q/a  t]*. 
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Figure  38.  Computed  steady-state  characteristic  of  the  inlet  orifice. 

Figures  39  and  40  show  the  steady-state  characteristics  for  the  vortex 
chamber  and  the  central  orifice  respectively.  The  average  tangential 
velocity  in  a vortex  is  approximately  proportional  to  the  square  root 
of  the  pressure  drop  across  the  vortex.  The  average  tangential  velocity 
in  a given  vovtex  chamber  increases  with  higher  flow  rates  and  lower 
velocities.  Therefore,  if  the  viscosity  decreases  (due  to  temperature 
rise)  the  flow  rate  must  decrease  to  maintain  a constant  average  tan- 


gential velocity  (or  equivalently,  the  pressure  drop), 
is  negative  for  the  vortex  chamber. 


Thus,  dQ/d T, 


AP=constant 


^Temperature  sensitivity  is  defined  for  a constant  pressure  drop. 


In  the  case  of  the  central  orifice  the  effect  of  inlet  swirl 
is  to  force  the  fluid  towards  the  cylindrical  wall,  thereby  reducing  cne 
effective  cross-sectional  area  of  che  orifice.  As  the  fluid  temperature 
increases  (accompanied  by  viscosity  decrease  and  swirl  increase)  the 
effective  orifice  area  decreases  and  the  resistance  increases;  the  result 
is  a negative  *Q/3t|aP  = constant  . 

At  low  pressure  drops  the  swirl  effect  is  small  in  both  the  vortex 
chamber  and  the  central  orifice.  In  this  case,  the  inlet  orifice  character- 
istic dominates  the  NSVR  resistance  characteristic,  resulting  in  an  overall 
positive  temperature  sensitivity.  At  higher  pressure  dropt  the  vortex 
chamber  and  the  central  orifice  characteristics  dominate,  resulting  in  a 
negative  temperature  sensitivity  for  the  NSVR.  If  the  inlet  orifice  can 
be  operated  in  the  pressure  drop  range  where  dC<j/dRe  is  negative  (possible 
only  for  L/lKl;  see  Figure  41),  the  inlet  orifice  characteristic  will  rein- 
force the  vortex  chamber  and  central  orifice  characteristics,  with  the  result 
that  the  overall  temperature  sensitivity  is  more  negative  than  with  the 
vortex  chamber  and  the  central  orifice  alone. 


Figure  41.  Discharge  coefficient  of  circular  short-tube  orifice 
with  sharp-edged  entry  (see  reference  18). 


One-Dimensional  Model  of  the  SVR 


4.3.1.  Model  Description 

The  SVR  geometry  is  conveniently  divided  into  five  regions  as 
shown  in  Figure  42  (see  Appendix  C for  mathematical  details).  Models  for 
regions  1,  2,  and  3 are  the  same  as  those  developed  for  the  NSVR.  A one- 
dimensional,  axisymmetric  flow  is  assumed  to  exist  in  region  4,  which  is 
contradictory  to  the  flow  field  observed  in  the  flow  visualization  study 
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(R*gions  2 & 4 


Figure  42.  Regions  of  SVR. 

described  earlier  (see  Section  4.1).  In  region  4 the  flow  field  is 
essentially  non-symmetric  about  a mid  plane  through  the  diffuser 
chamber.  But,  a non-uniform  velocity  profile  would  lead  to  partial 
differential  equations,  thereby  increasing  substantially  the  mathematical 
complexity  of  the  model.  The  simpler  one-dimensional  axisymmetric  assump- 
tion yields  a model  which  is  of  lesser  accuracy  but  greater  mathematical 
simplicity  than  the  partial  differential  equation  model. 

Regions  1 and  5 of  the  SVR  are  modeled  as  orifices  with  discharge 
coefficients  equal  to  that  of  a short- tube  orifice  with  a sharp-edged 
entry  (see  Figure  41). 


4.3.2.  Typical  Results 

Figures  43-46  show  comparisons  of  analytical  predictions  with 
typical  experimental  data  for  four  SVR  geometries.  Geometries  were  avail- 
able which  resulted  in  small  temperature  sensitivities  only,  although 
larger  temperature  sensitivities  are  possible  with  other  geometries.  The 
difference  between  the  experimental  and  predicted  values  of  the  flow  rate 
(at  constant  pressure  drop)  is  less  than  5 percent  of  reading  for  all 


but  the  low  flow  rate  range  in  all  four  cases  considered.  The  character- 
istics at  10*C  and  65.55*C  shown  in  Figures  43-46  indicate  the  temperature 
sensitivity  of  the  goemetries  over  a wide  temperature  range. 


Pressure  Drop 

Figure  43.  Comparison  of  predictions  of  one-dimensional  model 
with  experimental  data  of  SVR. 
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Figure  44.  Comparison  of  predictions  of  one-dimensional  model 
with  experimental  data  of  SVR. 
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Figure  45.  Comparison  of  predictions  of  one -dimensional  model 
with  experimental  data  of  SVR. 
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Figure  46.  Comparison  of  predictions  of  one -dimensional  model 
with  experimental  data  of  SVR. 


4.4. 


Two-Dimensional  Model  of  the  NSVR 


4.4.1. 


Model  Description 


The  steady-state  relations  used  in  the  one-dimensional  model  of 
the  NSVR  for  regions  1,  3,  and  4 are  also  used  for  the  two-dimensional  model. 

The  model  for  region  2 is  developed  by  adapting  the  momentum 
integral  approach  of  Wormley^  to  the  laminar  flow  case.  Appropriate 
non-uniform  velocity  profile  functions  are  assumed  for  the  radial  and 
tangential  velocities,  and  the  laminar  shear  stresses  (both  radial  and 
tangential)  are  computed  directly  from  the  velocity  profiles.  Integra- 
tion of  the  partial  differential  equations  along  the  vortex  chamber  height 
leads  to  a system  of  four  non-linear  ordinary  differential  equations  which 
are  solved  numerically.  The  model  includes  a "developing-f low"  region  as 
was  originally  employed  by  Wormley^.  In  the  developing  flow  region  the 
radial  velocity  in  the  core  is  non-zero  and  the  circulation  (i.e.,  rv^) 
is  constant.  When  the  flow  is  "developed,"  the  radial  velocity  in  the 
core  is  zero  and  the  circulation  is  not  constant.  Unlike  the  one-dimen- 
sional model,  the  shear  stresses  are  computed  directly  from  the  velocity 
profiles . 


The  differential  equations  are  locally  linearized  at  r = r^ 

(see  Appendix  D for  details)  to  eliminate  computational  difficulties  which 
result  from  the  singularity  at  the  outer  periphery  of  the  vortex  chamber 
(i.e.,  at  r=r-^).  An  approximate  solution,  valid  in  a close  neighborhood 
of  rj,  is  obtained  to  initiate  the  numerical  integration. 


4.4.2.  Typical  Results 

Figures  47-50  show  comparisons  of  model  predictions  and  typical 
experimental  data  for  different  NSVR  geometries.  The  difference  between 
experimental  and  predicted  values  of  flow  rate  (for  a constant  pressure 
drop)  is  less  than  15  percent  of  reading  in  all  cases. 

The  two-dimensional  model  predicts  the  temperature  sensitivity 
more  accurately  (see  Figures  33  and  48)  than  the  one-dimensional  model. 
The  improvement  in  the  accuracy  may  be  attributed  to  the  two-dimensional 
nature  of  the  velocity  profiles. 

Figures  51  and  52  show  the  radial  distribution  of  the  velocities 
and  the  boundary  layer  thickness  corresponding  to  the  NSVR  geometry  of 


^Wormley,  D.  N.  "Analytical  Model  for  the  Incompressible  Flow 
in  Short  Vortex  Chambers,"  Trans.  ASME,  J.  Basic  Engineering,  Series  D, 
Vol . 91,  No.  1 (June,  1969)  , pp.  264-276. 
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Figure  47  for  a flow  rate  of  0.6554x10'  mJ/s  (at  35  C) 
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Comparison  of  predictions  of  two-dimensional  model 
with  experimental  data  of  NSVR. 
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Figure  51.  Radial  distribution  of  the  radial  velocity  functions 
(u$ , us)  and  the  boundary  layer  thickness  (6)  as  predicted  by  the 

two-dimensional  model. 


Radial  Distance  (mm) 

Figure  52.  Radial  distributions  of  tangential  velocity  ana  pressute 
in  the  vortex  chamber  as  predicted  by  the  two-dimensional  model. 


The  radial  velocity  in  the  core  (i.e.  u^)  is  zero  at  about  r=5.08mm 

(see  Bigure  51),  indicating  that  all  the  flow  is  contained  in  the  boundary 
layer.  The  boundary  layer  thickness  (6)  increases  rapidly  in  the  developing 
flow  region.  The  thickness  (6)  increases  slowly  and  then  decreases  as 
the  flow  progresses  towards  the  center.  The  decrease  is  due  to  a rapid 
acceleration  of  the  fluid  particles  in  the  radial  direction  caused  by  the 
steep  pressure  gradient  existing  near  the  center  (see  Figure  52).  The 
tangential  velocity  (v§)  increases  linearly  in  the  developing  flow  region 
and  remains  relatively  constant  in  the  developed  flow  region. 


5.  DYNAMIC  BEHAVIOR  OF  NSVR 


This  section  summarizes  the  development  of  an  approximate  analytical 
model  which  describes  the  dynamic  behavior  of  the  NSVR  configuration,  and 
the  results  of  experimental  transient  response  studies  on  one  particular 
NSVR  geometry. 


lamic  Model  Development 


A detailed  discussion  of  the  dynamic  model  development  is  given 
in  Appendix  E.  The  dynamic  model  employs  the  one-dimensional  steady-state 
model  (Section  4)  with  appropriate  changes  made  to  include  primary  dynamical 
effects.  The  orifices  and  the  diffuser  chamber  are  modeled  as  purely  resis- 
tive elements.  A time  delay  T^  and  a first-order  lag  T2  characterize  the 
propagation  of  the  tangential  velocity  in  the  vortex  chamber.  The  resulting 
mathematical  model  is  nonlinear. 

Linearization  of  the  nonlinear  model  results  in  the  following 
transfer  function  relating  to  principal  variables  of  interest: 


APq  K3(1+T2D) 

= 1+T2EH-K4exp(-T1D) 

where  K3  and  K4  are  constants  defined  in  Appendix  E. 
used  for  small-signal  dynamic  analysis. 


Equation  2 can  be 


5.1.1. 


Determination  of  the  Time  Delaj 


and  the  Time  Constant  T-j 


The  time  delay  T^  and  the  time  constant  T2  are  estimated  by 
modeling  the  vortex  chamber  as  a "tube"  of  width  b and  length  L (see 
Figure  53).  The  time  delay  T^  is  assumed  to  be  the  time  taken  for  the 
tangential  velocity  information  to  travel  the  length  L.  That  is 
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Figure  53.  Tube-equivalent  of  vortex  chamber. 

where  v is  the  average  tangential  velocity  in  the  vortex  chamber  obtained 
from  the  steady-state  analysis. 

The  time  constant  T2  is  obtained  by  considering  the  inertance  I 
of  the  equivalent  tube  of  length  L and  the  overall  steady-state  resistance 
R of  the  vortex  resistor.  Therefore 


T2  = I/R 

The  inertance  I is  given  by  pL/bh. 


(4) 


5.1.2 


Computed  Step  Response 


The  dynamic  model  constants  (K3,K^,  Ti  and  T2)  are  computed  by 
using  average  flow  conditions.  Table  II  lists  numerical  values  for  an 
example  problem.  Flow  variables  (i.e.  velocities  and  pressure)  corres- 
ponding to  zero  flowrate  are  used  as  the  initial  conditions  for  the  dynamic 
simulation. 

Figure  54  shows  the  computed  transient  response  of  a NSVR  (see 
Table  II  for  geometry)  subjected  to  a step  input  of  137.834  kN/m^.  The 
overshoot  is  due  to  the  time  delay  T^  associated  with  the  tangential 
velocity  propagation.  Until  the  new  tangential  velocity  information 
reaches  the  central  orifice,  the  NSVR,  exhibits  a resistance  value  corres- 
ponding to  the  initial  operating  point.  For  a positive  step  input,  this 
delay  corresponds  to  a temporary  increase  in  the  resistance.  The  temporary 
decrease  in  flow  rate  caused  by  the  delay,  results  in  the  overshoot  observed 
in  Figure  54.  For  a negative  step  input,  there  would  exist  a corresponding 
undershoot. 
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Figure  54.  Computed  transient  response  of  the  NSVR  to  a step 


TABLE  II 

DYNAMIC  MODEL  PARAMETERS  - EXAMPLE 


Geometry  of  NSVR 

b = 3.810mm 

h = 8.890mm 

£ = 4.340mm 

r|  = 6.350mm 

T2  = 2.552mm 

j^2  = 10.160mm 

Model  Coefficients 


L = 32.367  mm 
V = 2.43  m/s 

I =793.02  kN/m  -sZ 
R = 3. 785x10"’  kN-s/m5 

Tj  = 13  milli-seconds 
T2  = 2 milli-seconds 
K3  - 0.8745  (Dimensionless) 

K/,  = 0.2796  (Dimensionless) 


5.2. 


Experimental  Study 


5.2.1.  Apparatus 

Figure  55  shows  a schematic  of  the  apparatus  used  in  the  experi- 
mental study  of  the  dynamic  behavior  of  vortex  resistors.  The  hydraulic 
system  used  was  the  same  as  that  used  for  the  steady-state  experiments. 

A quick  opening  solenoid  valve  [3/16  inch  (4.7625mm)  orifice  diameter]  was 
employed  to  generate  a transient  input  at  the  upstream  or  downstream  end 
of  the  vortex  resistor.  An  accumulator  and  a flow  bypass  path  (using  a 
needle  valve)  was  used  to  minimize  pressure  fluctuations  upstream  of  the 
solenoid  valve,  during  and  after  sudden  opening  of  the  valve.  A fixed 
resistor  (short-tube  orifice  L = 6.350mm,  D = 1.5875  mm),  was  connected  at 
the  downstream  end  of  the  vortex  resistor  to  simulate  a typical  load.  The 
instantaneous  pressure  drop  across  this  "load"  orifice  was  a measure  of  the 
instantaneous  flow  rate  at  the  exit  of  the  vortex  resistor.  Figure  56 
shows  the  dimensions  of  the  end  blocks  of  the  vortex  resistor  and  the 
pressure  transducer  locations. 
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Figure  55.  Schematic  of  transient  response  experiment, 
(solenoid  valve  was  located  at  the  downstream 
of  vortex  resistor  in  some  cases) 


5.2.2.  Instrumentation 

Kistler  Piezoelectric  Model  601A  transducers  with  Model  504a 
charge  amplifiers  were  used  to  measure  the  time  dependent  pressures  at 
the  inlet  and  exit  of  the  vortex  resistor.  The  long  time  constant  setting 
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Figure  56.  Vortex  resistor  block. 

on  the  charge  amplifiers  permitted  the  static  calibration  of  the  pressure 
measuring  system.  The  temperature  measurement  system  used  in  steady- 
state  experiment  was  used  to  measure  the  average  temperatures  during 
the  dynamic  response  tests.  The  temperature  of  the  hydraulic  oil  was 
relatively  constant  (within  + 0.5  C)  during  the  transient  response  tests. 

The  input  and  output  pressure  signals  (from  the  upstream  and 
downstream  of  the  vortex  resistor)  were  displayed  on  a Tektronix  Model 
502  dual-beam  oscilloscope.  Polaroid  pictures  were  taken  of  the  oscillo- 
scope traces . 


5.2.3. 


Experimental  Results 


Figure  57  shows  a typical  response  obtained  on  a NSVR  configura- 
tion. A positive  step  input  was  applied  by  suddenly  opening  the  solenoid 
valve  located  at  the  upstream  end  of  the  vortex  resistor  block.  The  rise 
time  of  the  pressure  "step"  input  was  of  the  order  of  10  milli-seconds ; 
this  was  the  shortest  rise  time  possible  with  the  available  components. 

No  discernable  pure  time  delay  or  first-order  lag  is  evident  on  the  tran- 
sient response  traces  of  Figure  57. 


1 div,  = 20  milli-seconds 

Figure  57.  Transient  response  of  a NSVR  to  a positive  step  input. 

Figure  58  shows  a transient  response  trace  taken  for  a some- 
what higher  pressure  step  input.  The  high  frequency  noise  in  the  input 
was  due  to  a short  term  oscillation  of  the  pressure  regulator. 

In  an  effort  to  sharpen  the  input  pressure  signal,  a negative 
step  was  used.  The  negative  step  input  was  generated  by  locating  the 
solenoid  valve  between  the  downstream  end  of  the  vortex  resistor  and  the 
loading  resistor.  The  rise  time  obtained  was  of  the  order  of  4 milli- 
seconds. Figures  59  and  60  show  transient  responses  for  two  different 
loading  resistors  at  the  downstream  of  the  vortex  resistor.  No  discernable 
dynamics  are  evident  from  these  traces.  The  resolution  obtainable  from 
the  scope  photographs  is  of  the  order  of  0.5  milli-second.  Therefore, 
delays  or  time  constants  smaller  than  0.5  milli-second  are  not  observable. 


div.  = 68.917  kN/n/ 


1 dlv.  = 20  milli-seconds 

Figure  58.  Transient  response  of  a NSVR  to  a positive  step  input 
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Figure  59.  Transient  response  of  a NSVR  to  a negative  step  input, 
(loading  orifice  diameter  = 3.175mm) 
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Figure  60.  Transient  response  of  a NSVR  to  a negative  step  input, 
(loading  orifice  diameter  = 1.5875mm) 


5.3.  Comparison  of  Theory  and  Experiment 

The  rise  time  of  the  experimental  input  was  considered  too  long  to 
permit  a direct  comparison  of  experimental  data  with  the  predicted  step 
response  of  Figure  54.  A terminated-ramp-type  input  signal  having  a 
rise  time  of  10  milliseconds  approximates  the  experimental  input  condi- 
tion. Figure  61  compares  the  experimental  and  predicted  dynamic  responses 
for  this  ramp  input.  The  model  predicts  the  basic  dynamic  behavior  well 
in  a qualitative  sense,  even  though  there  appears  to  be  some  error  in  T^ 
and  T2-  More  accurate  estimation  of  the  time  constants  T^  and  T2  may  be 
possible  using  a distributed  parameter  approach  similar  to  that  used  by 
Anderson-^,  but  the  procedure  presented  in  this  section  is  probably 
adequate  for  most  purposes. 


■^Anderson,  W.  W.  "A  Dynamic  Model  of  Vortex-Type  Fluid  Amplifier, 
(Ph.D.  dissertation,  Massachusetts  Institute  of  Technology,  1971.) 
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Figure  61.  Comparison  of  experimental  and  computed  transient  response 
of  a NSVR  (to  a terminated  ramp  input). 

6.  NOISE  CHARACTERISTIC 

6.1.  Measurement 

A schematic  of  the  experimental  set  up  of  the  noise  measurement 
system  is  shown  in  Figure  62.  The  accumulator  was  used  as  a fluid  power 
source.  This  was  done  to  minimize  the  noise  generated  by  the  hydraulic 

. j i a jj i j *.1.  > ..... 

source. 
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was  used  only  to  initially  charge  the  accumulator.  The  flow  rate  varied 
gradually  as  the  accumulator  was  discharged.  The  discharge  time  of  the 
accumulator  was  long  enough  to  permit  noise  measurements  at  various  flow 
rates. 

Kistler  Piezoelectric  transducers  with  associated  charge 
amplifiers  were  used  for  the  noise  measurements.  Polaroid  pictures  were 
taken  of  the  pressure  traces  on  the  oscilloscope. 


Discussion  of  Noise  Characteristic 


Figure  63  shows  the  noise  level  with  no  flow  through  the  resistor. 
This  noise  was  due  to  the  electronic  instrumentation,  namely  that  due  to  the 
charge  amplifiers.  i 


Upstream 


1 div.  = 0.0689  kN/ni 


Downstream 
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1 div.  = 2 milli-seconds 

Figure  63.  Noise  level  with  no  flow. 

Figures  64  show  the  noise  signal  (with  the  d.c.  part  filtered  out)  at 
various  flow  rates.  Three  different  geometries  were  tried  and  each 
'esulted  in  similar  noise  traces.  The  vortex-resistor-generated  noise 
appears  to  be  well  below  that  of  the  electronic  instrumentation.  The 
instrumentation  noise  level  corresponded  to  a signal-to-noise  ratio  of 
about  900.  A spectrum  analysis  of  the  ncise  signal  was  considered 
unnecessary  in  view  of  the  fact  that  the  resistor-based  noise  level  was 
so  low. 


Upstream 


1'  div.  = 0.0689  kN/m 


Downstream 


1 div.  = 2 milli-seconas 
(a)  Flow  Rate  = 2.29xl0-6  m /s 
Figure  64.  Noise  level  of  NSVR  at  various  flow  rates, 
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(b)  Flow  Rate  = 5 06x10"  m /s 
Figure  64.  Noise  level  of  NSVR  at  various  flow  rates 
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(c)  Flow  Rate  = 10.91x10"^  m^/s 
Figure  64.  Noise  level  of  NSVR  at  various  flow  rates. 

7.  CONCLUSIONS  AND  RECOMMENDATIONS 

7.1.  Conclusions 

This  study  has  established  that  (1)  temperature-insensitive 
vortex  resistors  can  be  built  in  practical  sizes,  and  (2)  temperature 
compensation  for  fluid  resistors  of  the  orifice  and  capillary  tube  types 
can  be  provided  by  means  of  vortex  resistors  and  (3)  vortex  resistors 
(operating  on  hydraulic  oil)  behave  primarily  as  resistors  with  no 
appreciable  reactive  dynamics. 

The  empirical  models  discussed  in  Appendix  A can  be  employed 
for  designing  vortex  resistors  having  a desired  temperature  sensitivity. 

The  accuracy  of  these  models  cannot  be  insured  for  operation  outside  the 
range  of  the  data  base. 

The  one-dimensional  analytical  models  developed  for  the  NSVR  and 
SVR  configurations  can  also  be  used  for  design  purposes.  These  models 
are  expected  to  be  more  accurate  over  a wider  range  of  geometries  and  flow 
conditions  than  the  empirical  models.  The  empirical  models  do  not  include 
all  the  important  geometrical  parameters  and  are  limited  by  constraints 
placed  on  the  establishment  of  a data  base  through  experimentation.  The 
analytical  models,  include  most  of  the  variables  believed  to  be  of  signifi- 
cance. 
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The  two-dimensional  model  developed  for  the  NSVR  is  more 
accurate  than  the  one-dimensional  model  for  some  geometries,  but  it  is 
not  as  efficient  computationally  due  to  the  numerical  integration  required. 

The  dynamic  model  developed  for  the  NSVR  is  only  an  approximate 
dynamic  characterization  of  the  resistor,  but  it  is  sufficiently  accurate 
for  most  system  dynamics  needs.  Model  predictions  as  validated  by  trans- 
ient response  requirements,  suggest  that  the  dynamics  of  a vortex  resistor 
are  relatively  unimportant  compared  to  the  dynamics  of  other  principal 
elements  in  a typical  fluid  system  which  incorporates  the  resistor. 

Limited  noise  studies  have  shown  that  a signal-to-noise  ratio 
of  900  or  more  is  possible  with  a vortex  resistor  operating  in  the  laminar 
flow  regime. 


7.2.  Recommendations 


It  is  recommended  that  HDL  give  consideration  to  supporting 
continuing  studies  as  discussed  below; 

(1)  Steady-State  Empirical  Models;  The  empirical  models 
developed  for  the  SVR  and  NSVR  (low  resistance  flow  direction)  are 
limited  in  usefulness  to  a relatively  small  range  of  geometries,  tem- 
peratures, and  flow  conditions.  The  most  serious  limitation  of  the 
data  base  used  to  develop  the  empirical  models  is  the  small  temperature 
range  covered.  The  lack  of  a cooler  on  the  hydraulic  test  stand  prevented 
measurements  below  35  C.  The  kinematic  viscosity  of  MIL-5606A  hydraulic 
oil  varies  by  a factor  of  2 for  the  temperature  range  used  (35  to  57.2°C). 
Military  specifications  usually  call  for  an  operating  temperature  range 

of  -45.55  to  65.55  C,  or  a kinematic  viscosity  variation  of  500  to  8 
centistokes  for  MIL-5606A  hydraulic  oil.  Additional  data  which  cover 
this  wider  temperature  range  are  needed  to  validate  the  analytical  models 
and  to  improve  the  empirical  models.  A cooler  could  be  installed  on  the 
OSU  test  stand  at  relatively  low  cost. 

(2)  Steady-State  Analytical  Models;  Velocity  profiles  in  the 
vortex  chamber  should  be  measured  to  validate  more  fully  the  flow  field 
assumptions  which  under ly  the  two-dimensional  model. 

An  improved  diffuser  model  should  be  sought  for  the  SVR 
configuration.  The  diffuser  model  for  the  SVR  should  take  into  account 
the  secondary  reverse  flow  in  the  diffuser.  A momentum  integral  approach 
can  be  used  for  this  purpose. 

(3)  Applications  of  the  Analytical  Models:  Algorithms  and 

computer  programs  could  be  developed  to  make  use  of  the  analytical  models 
as  design  tools.  Optimization  capability  can  be  incorporated  into  these 
design  aids  with  relative  ease. 
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(4)  Dynamic  Response  Characteristics:  The  dynamic  response 

experiments  conducted  to  date  have  established  that  the  time  constants 
associated  with  a typical  vortex  resistor  (vortex  chamber  diameter  = 12.7mm) 
are  less  than  0.5  milliseconds.  More  precise  determination  of  these  time 
constants  would  be  difficult  using  typical  size  geometries.  A more  accu- 
rate dynamic  characterization  could  be  made  by  time-scaling  transient 
response  measurements  on  a large  scale  model. 

A distributed  parameter  model  similar  to  Anderson's  turbulent 
flow  model-^  could  be  developed  to  more  accurately  predict  the  vortex 
chamber  dynamic  characteristics  for  the  laminar  flow  case.  This  model 
would  yield  values  for  the  time  constants  T^  and  T2  to  be  used  in  the 
lumped-parameter  model. 

(5)  Noise  Characteristics:  This  study  has  established  only  a 

lower  bound  (of  900)  for  the  signal-to-noise  ratio  of  a vortex  resistor. 
Improved  measurement  techniques  and  alterations  in  the  test  loop  would 
allow  more  precise  signal-to-noise  measurements  if  such  information  is 
needed . 
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APPENDIX  A. 


STEADY-STATE  EMPIRICAL  MODELS 


This  appendix  describes  the  development  of  empirical  models  of 
symmetric  (SVR)  and  non-symmetric  (NSVR)  vortex  resistors. 


A.l. 


Basic  Model 


The  geometries  of  the  NSVR  and  the  SVR  are  described  by  six 


b,  h,  r^,  r2>  and  l Due  t0  symmetry. 


length  parameters  (see  Figure  4); 

the  diffuser  chamber  dimensions  are  the  same  as  those  of  vortex  chamber  for 
a SVR.  The  steady-state  characteristic  of  a fluid  resistor  may  be  described 
by  a functional  relationship  between  the  volumetric  flow  rate  (for  incom- 
pressible flow  situations)  and  the  pressure  drop  across  the  resistor.  This 
functional  dependence  involves  the  geometry  of  the  resistor  and  relevant 
physical  properties  (density  and  viscosity  in  this  case)  of  the  fluid. 

In  the  case  of  the  vortex  resistor  this  relationship  may  be  expressed  as 
follows : 


Q - Q (b , h,  ^ 1 9 ^1*  r2  ’ ^2’  * P*  U» ) 


(A.l) 


For  purpose  of  scaling,  it  is  necessary  to  normalize  the  variables.  The 
following  dimensionless  parameters  may  be  defined  based  on  the  variables 
of  equation  (A.l): 
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The  choice  of  r as  the  characteristic  length  is  arbitrary.  Equation 
(A.l)  may  be  replaced  by 
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The  function  f can  be  any  function  that  fits  the  experimental  data  within 
the  desired  accuracy.  However,  the  choice  of  the  function  f is  governed 
by  mathematical  simplicity  and  ease  of  numerical  manipulation.  From  this 
point  of  view,  a reasonable  functional  form  would  be  as  follows: 
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where  a^,  a^,  , a^  are  constants.  These  constants  can  be  determined 

by  linear  regression  using  experimental  data. 


A. 2.  Data  Base 

Ideally  the  experimental  data  to  be  used  for  an  empirical  model 
development  should  envelop  with  a fine  grid  the  entire  range  of  useful 
geometries  and  flow  situations.  However,  practical  considerations  make 
such  an  exhaustive  data  acquisition  infeasible.  Therefore,  a limited 
number  of  geometries  were  selected  and  tested  under  a few  judiciously 
selected  flow  situations.  Various  factors  were  involved  in  the  selection 
of  the  geometries  and  the  operating  conditions.  These  factors  are  discussed 
below. 

For  the  NSVR  configuration,  24  geometries  were  selected  close  to 
the  geometrical  proportions  of  the  Lee  Visco  Jet  #62AVL1  (see  reference 
?)  which  was  observed  to  have  negative  dR/^T  while  operating  on  MIL-5606A 
hydraulic  oil.  Of  the  five  geometrical  ratios  (i.e.,  b/rj,  h/r.,  £,/r. , ^/r ^ , 
and  ?2/ri),  only  three  were  considered  as  having  a strong  inf luenceonthe 
vortex  characteristics.  These  ratios  were  b/r^,  h/r^  and  r2/r^.  This 
simplification  was  necessary  to  limit  the  number  of  geometries  to  a manage- 
able number. 

The  vortex  chamber  radius  r^  was  selected  to  be  0.35mm.  Although 
smaller  values  of  r would  likely  be  used  in  most  fluidic  applications, 
the  dimensional  accuracy  afforded  by  the  available  machining  equipment 
prohibited  a smaller  value  of  r^  for  the  experiments  reported  here.  Figure 
A.l  shows  the  grid  locations  of  the  24  geometries  used.  The  dimensions  of 
these  geometries  are  given  in  Table  A.l.  In  each  of  the  six  vertical  planes 
of  Figure  A. 1 four  geometries  were  selected,  thereby  allowing  uniform  spanning 
of  the  geometrical  space.  The  geometries  were  staggered  to  decouple  the  space 
variables.  Otherwise,  the  selection  of  the  geometries  was  arbitrary.  The 
data  obtained  on  24  geometries  for  a NSVR  operating  in  the  high  resistance 
direction  of  flow  is  considered  adequate  for  development  of  the  NSVR  model. 

Additional  experimental  data  were  obtained  for  the  SVR  configura- 
tion and  for  flow  in  the  low  resistance  direction  of  the  NSVR.  The  ratio 
r2 /r i vas  held  constant  for  these  data.  These  data  were  primarily  intended 
for  evaluation  of  the  analytical  model.  Therefore,  only  a few  geometries 
were  selected  from  the  available  laminations  (see  Table  A. II.).  The  number  of 


geometries  used,  i.e.  8,  is  not  sufficient  to  develop  an  empirical  model 
valid  over  a wide  range  of  operating  conditions.  However,  the  empirical 
models  developed  using  these  limited  experimental  data  can  provide  at 
least  approximate  estimates  of  the  steady-state  characteristics  of  the  SVR 
and  of  the  NSVR  operating  in  the  low  resistance  flow  direction. 


2Technical  Hydraulic  Handbook,  The  Lee  Company,  Westbrook, 
Connecticut,  1$71. 
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TABLE  A. II.  GEOMETRIES  OF  THE  VORTEX  RESISTORS 
(CONSTRUCTED  BY  AIR  RESEARCH  CORPORATION) 


Geometry 

Number 

b 

u m 

h 

mm 

£1 

mm 

dl 

mm 

d2 

mm 

*2 

mm 

1 

0.7213 

1.524 

7.5660 

6.9240 

3.048 

0.254 

2 

0.7213 

1.524 

7.5660 

6.9240 

3.048 

1.778 

3 

1.2115 

0.762 

5.7658 

7.1988 

3.048 

0.254 

4 

1.2115 

0.762 

5.7658 

7.1988 

3.048 

1.778 

5 

0.7213 

0.762 

7.5660 

6.9240 

3.048 

1.778 

6 

0.7213 

1.524 

7.2771 

7.1988 

3.048 

0.254 

7 

0.9321 

1.524 

7.2771 

7.1988 

3.048 

0.254 

8 

0.9321 

1.524 

7.2771 

7.1988 

3.048 

1.778 

t-  l 


S 


A. 3.  Results 


Table  A. Ill  lists  the  constants  (a^s)  for  the  NSVR  and  the  SVR 
for  various  flow  situations.  Figures  A.2-A.6  show  correlations  of  the 
experimental  data  with  the  empirical  models.  The  empirical  predictions 
are  within  ten  percent  for  most  of  the  cases.  In  cases  2-4  the  maximum 
error  between  model  and  data  is  less  than  20  percent  and  in  case  1 the 
error  is  within  30  percent. 


From  equation  (A. 3),  with  ag  - ay  = 0, 


0 = 


rx4  + a5  AP 


1 

2+a5 


L-i  <!/2  <;/3 


2-fflc 

v 5 


(A.4) 


Differentiating  with  respect  to  the  kinematic  viscosity 


*ai 

dv  | 

AP  = constant 


(A. 5) 
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where 


1 


1 


r^S 


2+a, 


AP 


(7-)  2 (7-)  3 (7s) 
ri  ri  1 


*2/4 


From  Equation  (A. 5)  and  Table  A. Ill  it  can  be  observed  that  the  sign  of 
the  viscosity  sensitivities  (dQ/dv)  for  the  NSVR  (in  reverse  flow)  and 
the  SVR,  are  negative  for  low  Reynolds'  number  and  positive  for  high 
Reynolds'  numbers.  The  sensitivities  are  always  negative  for  the  forward 
flow  direction  of  the  NSVR.  A positive  value  of  Aq/Sv  implies  that  the 
flow  rate  increases  with  viscosity.  A small  magnitude  of  AQ/dv,  as  in 
the  case  of  the  NSVR  (high  resistance  direction  of  flow)  for  the  *2^rl  = 
0.402  geometries  indicates  that  the  resistor  is  nearly  temperature 
insensitive. 


xioa 


Figure  A. 2.  Empirical  model  of  NSVR  operating  in 
high  resistance  flow  direction. 
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APPENDIX  B. 


ONE-DIMENSIONAL  MODEL  OF  NSVR 


A one-dimensional  model  is  derived  in  this  appendix.  The 
modeling  of  the  NSVR  involves  division  of  the  resistor  space. 

Figure  B.l  shows  the  four  regions  of  the  NSVR  which  can  he 
individually  modeled.  The  equations  derived  here  are  similar  to  those 
of  reference1^  but  with  modif ications  for  operation  in  laminar  regime. 
Also  reference^  does  not  include  region  4 of  Figure  B.l. 


Figure  B.l.  Regions  of  NSVR 


12 

Bichara,  R.  T.,  and  P.  A.  Orner.  "Analysis  and  Modeling  of 
the  Vortex  Amplifier,"  Trans . ASME,  J- Basic  Engineering,  Series  D,  Vol . 91, 
No.  1 (December,  1969),  pp.  755-763. 
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Further,  due  to  operation  in  laminar  regime  the  discharge  coefficient  of 
inlet  orifice  is  not  treated  as  constant  but  is  allowed  to  vary  with  the 
Reynolds  number.  The  basic  assumptions  are: 

1.  The  fluid  is  incompressible  and  Newtonian, 

2.  The  flow  is  laminar  throughout  the  device, 

3.  No  heat  transfer, 

4.  All  surfaces  are  smooth  and  geometrically  exact, 

5.  No  body  forces, 

6.  The  operating  fluid  properties  are  assumed  constant  through- 
out the  flow  field. 

Additional  assumptions  are  described  in  the  sections  which  follow. 


B.l.  Region  1 

Figure  B.2  shows  region  1 of  the  NSVR  (same  as  region  1 of  SVR) . 


Figure  B.2.  Region  1 of  vortex  resistors.  (SVR  and  NSVR) 


B.1.1.  Assumptions 

1.  A shear  stress  based  on  the  flat  plate  Blassius  profile  is 
assumed  to  exist  on  the  annular  surface  of  region  1. 

2.  The  radial  and  tangential  velocities  and  pressure  are  uniform 
over  the  annular  surface  of  region  1. 

3.  The  discharge  coefficient  of  the  inlet  passage  is  the  same 
as  that  of  an  equivalent  circular  short  tube  orifice  with  sharp-edged 
entry  even  though  the  downstream  conditions  are  not  the  same  for  the  two 
cases. 

4.  The  width  of  the  inlet  orifice,  b,  is  much  smaller  than  the 
diameter  of  vortex  chamber  (i.e.,  b«2r^). 

B.1.2.  Governing  Equations 


Flow  Rate: 


Q = Cdlbh“\y2(pt  ^ pt)/p 


(B.l) 
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Conservation  of  i pr^  2 

angular  momentum  bhcj [ Q*  ‘ Pr]Qvl  = ^rrri  ^Tti  (B.3) 

Shear  Stress: 

Ttl  = 0.664  ^v11,5  / Vvxml  (B.4) 

where 

-i  /rrb 

- COS  I 

\ rl 

Orifice  Discharge  Coefficient: 

5 6 ,.,/M,., 

Cdl  '£l  £l  Cl3  <l08  Re‘‘e)  \~J  <B-5) 


where 


40  . . 2bh 

Rede  = irdv  and  de  = b^h 
e 

C(1,1)=0. 209285657990720  D 00  ,C(2,1)=0. 791788479941879  D-02 
C(3,l)=0. 226876244315706  D 00  ,C(4, 1)=-0.507148451077533d-01 
C(5,l)=-. 249327806165838  D-02  ,C(1,2)=-0.795026948224413D  00 
C(2,2)=0. 142923814768118  D 01  ,C(3,2)=-0. 10706294747 1109D  01 
C(4,2)=0. 327389281735691  D 00  ,C(5 ,2)=-0. 322173197963765D-01 
C(l,3)=0. 632542509986194  D 00  ,C(2,3)=-0. 127600317094214D  01 
C(3,3)=0. 934492443943111  D 00  ,C(4,3)=-0.286766619747770D  00 
C(5,3)=0. 304900339103475  D-01  ,C(l,4)=-0. 15149279067314 ID  00 
C(2,4)=0. 312611655321839  D 00  ,C(3,4)=-0.229092531851455D  00 
C(4,4)=0. 707338263930315  D-01  ,C(5,4)=-0.766974040642524d-02 
C(l,5)=0. 148165042552580  D-01  ,C(2,5)=-0.307706419128231D-01 
C(3,5)=0. 225854454114608  D-01  ,C(4,5)=-0.700958813169788D-02 
C(5,5)=0. 768051230807445  D-03  ,C(l,6)=-0.518515906302915D-03 
C(2,6)=0. 107854148254401  D-02  ,C(3,6)=-0.792969865409395D-03 
C(4,6)=0. 247201065756165  D-03  ,C(5,6)=-0.272806432199908D-04 
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i 


Equation  (B.5)  is  an  empirical  model  developed  from  the  experimental  data 
of  reference^®  and  is  accurate  to  within  5 percent  error  in  the  range 

10  <Rede  2=2200 


0.5  < * /d  ^12. 
1 e 


Figure  B.3.  Region  2 of  vortex  resistors.  (NSVR  and  SVR) 


B.2. 

B.2.1. 


B.2.2. 


Region  2 
Assumptions 

1.  Flow  is  symmetric,  see  Figure  B.3. 

2.  Flow  is  one-dimensional,  i.e.,  pressure  and  velocities 
(radial  and  tangential)  are  functions  of  only  the  radial 
distance. 

3.  Friction  in  the  radial  direction  is  neglected. 


Governing  Equations 


■ g ^ — — 

Radial  Momentum: 

du  v2 
u dr  ' r 

1 dj) 
p dr 

(B  .6) 

Tangential  Momentum: 

dv  , uv 
u — + — 
dr  r 

1 9rt2 
~ p Ay 

(B.7) 

Continuity: 

h (ru) 

= 0 

(B.8) 

Shear  Stress: 

't2Uf 

h 

(B.9) 

l®Reid,  K.  N.  "Static  Characteristics  of  Fluid  Amplifiers," 
Proceedings  of  the  Fluid  Power  Research  Conference,  Oklahoma  State  Univer- 
sity, July  25-26,  1967. 
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This  shear  stress  value  is  twice  that  for  fully  developed  laminar  flow 
between  parallel  plates. 


B.3.  Region  3 

Figure  B.4  shows  region  3 of  NSVR. 


B.3.1.  Assumptions 


1.  Radial  velocity  has  a discontinuity  at  region  3's  boundaries 
in  vortex  and  diffuser  chambers.  The  radial  velocity  is  zero  in  region  3, 

2.  Flow  is  axisymmetric , 

3.  The  rotary  motion  of  fluid  is  that  of  a rigid  body, 

4.  The  tangential  shear  stress  on  the  cylindrical  surface  is 
given  by  ths  flat  plate  boundary  layer  stress. 


B.3.2.  Governing  Equations 


Tangential 

velocity 

V = — V2 
r2  2 

(B. 10) 

Pressure  j 

distribution 

dp  2 

_ _ Pv 
dr  ~ r 

(B. 11) 

Axial  velocity 

w =-y/2(pu  - p3/p 

(B. 12) 

Boundary  condition-)  „ . 1 2 

« r - r,  < Pu  ' P2  + 2 PU2 

(B. 13) 
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Flow  rate 


angular  momentum 


^ ^d2  f 2nr  w dr 

0 

3 pr2Qv2  - pr2Qv3  = l^ 


(B.14) 


(B. 15) 


where 


t , = 0.664p, 


Region  4 


1 VVXm2  > xm2  ' 


B.4.1.  Assumptions 

1.  The  diffuser  chamber  is  assumed  to  be  large  enough  to  result 
in  small  velocity  heads. 

2.  The  diffuser  efficiency  is  assumed  to  be  0.5.  That  is  50 
percent  of  the  kinetic  energy  is  recovered  in  the  diffuser.  This  figure 
(50  percent)  is  selected  on  the  basis  of  the  axial  momentum  equation 
(written  for  the  control  volume  comprising  of  region  4 of  NSVR)  assuming 
potential  vortex  and  neglecting  axial  velocity  terms.  The  resulting 
axial  momentum  equation  is 


1 2 r *2  a2  , o 

-2>v3  L1  ’ K,  ’ ^ M 


which  gives  a pressure  recovery  of  approximately  50  percent  for  an  area 
ratio  (A  /A  ) of  0.2  (used  in  the  experiments).  (Assumptions  of  rigid  body 
rotationancl  constant  tangential  velocity  lead  to  unbounded  pressure 
recovery  for  increasing  diffuser  dimensions.) 


B.4.2.  Governing  Equations 


p0  ■ "j  + 4 “ 


2 4 + V3 

it  r 

2 


(B. 16) 


Equation  (B.l)  through  (B.16)  describe  the  mathematical  model  of 
the  NSVR  in  laminar  regime  of  operation.  These  equations  may  be  rearranged 
for  numerical  solution  as  follows: 
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From  equation  (B.l): 


Pi  = Pi  - 


JL 


Q2. 


2(Cdlbh)2 

Combining  equations  (B.3)  and  (B.4): 
1.5 


Vl 


+ a2v1  - a3  = 0 


where 


aL  ==  1.328nr1  hp 

a2  = P^Q 

2 

a3  s pr^Q  /bhCdl. 
From  Equation  (B.8): 

u2  = rlul/r2 • 

From  Equations  (B.7)  - (B.9): 


V2  = 


rl/r2 

exp  & h$I 


where 

2 

X s 24  vr1/h1 |u^ | . 
From  Equations  (B.6)  - (B.9): 


puL 

p2  = P1  2 \ ~^~2 


-1 


PVj^ 

exp (X ) 


oo  i f . . 2i-2-> 

X { l-(r2/r1)  } 

i!  (2i-2) 


(B.17) 


(B. 18) 


(B. 19) 


(B.20) 


(B  21) 
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Combining  Equations  (B.10)  - (B.14)  gives: 


q-5  ^ 


!nr22  '2 

2 L 


2^p2~p3^  + 


f (P2-P3)+u22 


(B . 22) 


From  Equation  (B.15): 


V-l-V-j  - 


0.664nvr2^2 

Q\Z2vxm2 


(v2  + v3> 


(B. 23) 


Numerical  Computations  Algorithm 


Given  the  geometry,  fluid  properties  and  flow  rate  Q the  pressure 
drop  &P=Pi  - pQ  may  be  computed  as  follows: 

Steps:  1.  Calculate  u^  from  Equation  (B.2),  and  p^  from 

Equation  (B.17)  by  assuming  any  arbitrarily  chosen 
value  of  p. . 

2.  Solve  Equation  (B.18)  for  v..  Due  to  nonlinearity  it 
will  involve  a numerical  search. 

3.  Calculate  v2  from  Equation  (B.20),  and  p2  from 
Equation  (B.21).  Note  that  the  series  in  Equation 
(B.21)  converges  to  a desired  accuracy  within  a finite 
number  of  terms  (less  than  50  terms). 

4.  Numerically  solve  Equation  (B.22)  for  p,,  and  Equation 
(B.23)  for  V3. 

5.  Compute  p from  Equation  (B.16)  and  calculate 


APPENDIX  C.  ONE-DIMENSIONAL  MODEL  OF  SVR 


Figure  C.l  shows  the  five  regions  of  SVR.  The  models  or  regions 
1,  2,  and  3 of  SVR  and  NSVR  (described  in  Appendix  B)  are  identical.  Only 
the  models  of  regions  4 and  5 are  described  below. 


C.l.  Region  4 

Figure  C.2  shows  region  4 of  SVR. 


Figure  C.2.  Region  4 of  SVR. 


1.  Flow  is  assumed  one -dimensional  and  axisymmetric . 

2.  A fraction  of  the  angular  momentum  at  the  inlet  of  the  diffuser 
chamber  is  recovered  as  potential  energy.  The  remainder  is  dissipated  com- 
pletely in  the  diffuser  chamber. 

3.  The  shear  stress  in  the  radial  direction  is  neglected.  Blassius 
solution  of  flat  plate  boundary  layer  is  used  for  the  tangential  shear  stress 
used  in  the  model. 


C.1.2.  Governing  Equations 

Radial  Momentum: 

Tangential  Momentum: 

Continuity: 

Shear  Stress: 

Boundary  Conditions : 

C.2.  Region  5 

Figure  C.3  shows  region  5 of  SVR. 


du  v^  _ _1  d£ 

dr  r p dr 


dv  uv  1 dT 

^dr  r p 3y 


dr 


(ru)  = 0 


+ h = 
y=~  2 


12pi 


+ i P (v. 


'3  =P3^2PVV3  *V3 


’2) 


v(r1)  = 0 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.3) 

(C.6) 
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C.2.1.  Assumptions 

1.  The  tangential  momentum  at  the  annular  inlet  of  region  5 is 
assumed  to  be  dissipated  completely  due  to  friction  at  the  annular  wall. 

2.  An  orifice  characteristic  may  be  used  for  the  exit  channel. 


C.2.2.  Governing  Equations 

Q = Cd3bh  'n/2(P4"P0)/P  (C.7) 

where  C^-j  = C^  and  is  calculated  from  Equation  (B.5). 


Equations  (B.l)  through  (B.15)  and  (C.l)  through  (C.7)  describe 
the  complete  SVR  model.  Equations  (B.17)  through  (B.23)  are  valid  for 
this  case.  Additional  equations  needed  for  calculations  of  p^  and  pQ  are 
developed  below:  From  Equations  (C.2)  - (C.4): 


r2/r 

exp 

^(r2  - r221 

2 2 

L rl  J 

(C.8) 


From  Equation  (C.8)  it  can  be  observed  that  the  boundary  condition  (C.6) 
cannot  be  satisfied  exactly.  Therefore  instead  of  requiring  v(r^  to  be 
zero,  it  was  allowed  to  be  small  enough  so  that  the  associated  angular 
momentum  was  a negligible  fraction  of  the  total  momentum.  Therefore 


where  v(r^)  = 0.000254  m/s  (assumed). 

Combining  Equations  (C.l),  (C.3),  and  (C.8)  gives  upon 
simplification: 


(C.9) 


1 

2 


OO 

T 

i=2 


f-x)1  fi-(Vri)21-2) 


i ’ (2i-2) 


(C.10) 
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From  Equation  (C.5) 


Po  = P4 


2(Cd3bh)‘ 


(C.ll) 


Equations  (B.17)  - (B.23)  and  (C.5),  (C.9),  (C.10),  and  (C.ll) 
constitute  the  SVR  model  for  numerical  calculations.  These  equations  are 
solved  sequentially  to  determine  the  pressure  drop  across  the  resistor  for 
a given  flow  rate. 


APPENDIX  D.  TWO-DIMENSIONAL  MODEL  FOR  NSVR 


In  Appendix  B a steady-state  one-dimensional  model  was  derived 
for  the  NSVR.  The  velocity  profiles  in  region  2 (i.e.,  the  vortex  chamber) 
were  assumed  to  be  uniform  across  the  chamber  height.  In  this  appendix 
this  restriction  has  been  relaxed  and  the  velocity  profiles  are  allowed 
to  be  nonuniform.  Appropriate  equations  are  derived  for  region  2 based 
on  Von  Raman' s momentum  integral  technique.  Equations  derived  earlier 
in  Appendix  B for  regions  1,  3,  and  4 remain  the  same. 


D. 1 . Region  2 

Figure  D.l  shows  region  2 of  NSVR. 


Figure  D.l.  Region  2 of  NSVR. 


D.1.1.  Assumptions 

1.  Vortex  chamber  is  thin,  i.e.,  h/r^  «1. 

2.  Flow  is  axisymmetric . 

3.  Flow  is  symmetric  about  xx  (see  Figure  D.l). 

4.  The  radial  velocity  in  the  core  is  zero  when  the  flow  is 
fully  developed. 


D.1.2.  Governing  Equations 


Continuity: 


ftu  u ?Vw 

?)r  r *y 


(D.l) 


Boundary  Layer 
Momentums : 


2 

du  v | 3u 
3r  ' r ^y 


1 dp  ^ 

r*-  + v 

P dr 


(D.2) 
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Core  Momentum: 


r 


(D.4) 


u 


Sv 


Boundary  Conditions: 


Velocity  Profiles: 


J6  dr 
dv6 


p dr 


^6 


U6V6 


= 0 


dr  r 

At  y = 0,  u = v = w = 0 
At  y = h/2,  w = du/dy  = dv/dy  = 0 
At  r = rp 


(D.5) 


u = U]_ 

V = 

w = 0 
P = Pi 


u = Ug  ( r ) f (T\)  + us(r)  g(Tl) 
v = v6  (r)f  (T|)  J 

u = ufi  (r) 1 

\ 6 < y < h/2 

v = v6  (r)  I 

H s y/6  J 


(D.6) 

0 < y < 6 

(D.7) 


Integrating  Equations  (D.l)  - (D.3)  using  the  velocity  profiles 
of  (D.7)  and  the  boundary  conditions  of  (D.6)  gives  the  following  system 
of  ordinary  differential  equations  (note  equations  D.4  and  D.5  are  also 
included) : 


[(35-l)u6-tar4usl  £ + £*5-1)6+!]  -1  + [«46]  =’7  [(“s'1) 


du 


6u6  + 1 u6  + *4  6 us] 


(D.8) 


[<«r«5>  u52  + (2  of2  - 0f4)  uAus  +a3ug2]  ^ + [(^l-V  6u6  + 

du  _du  . . c 

VUs]  —■  + [(^2^4>  6u6  + to36u.]  d7  + [p3  df  " 7 


I5lv62  ‘ “6  2 * (2*2^V  U6US  ' Vs2]  ‘ T (fo  u6 
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dr  + L(“l  6u6  +“26us>/v6J  dT  + 

[(«2  ' «4>3  dT  = I E(23,1^5>- 

u6  - (2“2  - «4>  u3  * I fo 


du6 

dr  + 

IpJdr 

dv*  = 
dr 

-u6v6/r 

(D. 10) 

(0.11) 

♦ 

(D . 12 ) 


where 


r1  2 r1 

= / f^dT)  , a 2 = / fgdTl 

o o 


1 2 1 

a3  = f « d1[l  » a4  = / gdi:i 
o o 


“s  - / £<r"  • £o'  - 3r 


T1  = 0 


' _ is 

3o  dTl 


T)  = 0 


Profile  functions  f and  g may  be  selected  arbitrarily  to  satisfy  the 
boundary  conditions  [f(0)=0,  f(l)=l,  f (0)>0,  f'  (1)=0  and  g (0)=g(l)=0 
g'(0)>0,  g'(l)=0].  The  following  functions  satisfy  these  conditions; 


fOD  = c[j  - {Hc/(c_1)}] 

gOt)  = b pKl-H)  H 


CD. 13) 

4 

(0.14) 


There  exists  a singularity  at  r=rlf  for  the  system  of  differential 
equations  (D.8)  - (D.12).  Approximate  solutions  for  6,  us  and  u$  may  be 
obtained  by  locally  linearizing  the  differential  equations.  The  resulting 
solution  is  valid  in  the  neighborhood  of  the  singularity. 


D.2.  Approximate  Solution  Near  the  Singularity 

Forcing  u§  to  zero  Equations  (D.9)  -(D.ll)  give  respectively 


a 


2 d6  ^ 0 c dus  6 dp 

-u  — + uo  -T-  + — -r*- 

3 s dr  J s dr  o dr 


s dr  p dr 


6 r 

= rL  1 


d6 


v6  "a3Us  J - 6 go'  u£ 


du 


(*2^V  us  d7  + (*2^4)6  dT 

(2v2^4>  vf  ' 


0 !' 


6u  - + 6us 

s 6 r 8 


1 dp  6 

p dr  r 

Assuming  potential  core  flow.  Equation  (D.17)  gives 

P(r)  = Pi  + f <rlvl)2  " ~2> 


(D.15) 


(D. 16) 
(D.17) 

(D. 18) 


Rearranging  Equation  (D.16) 


<“2^V  dr  (6U8>  = ‘ r 


2a2^4 


6u 


vf  1 0L2 

_ 6us  (D.19) 


At  r=rlf  ug  = 0 but  d(6us)/dr  is  not  negligible.  Therefore, 
neglecting  6ug  terms  in  Equation  (D.19)  gives 


d vf  ' 

^2-a4>  d?  (6us>  = • ~T~ 


(D. 20) 


Neglecting  the  radial  friction  in  Equation  (D.15)  and  combining  with 
Equation  (D.18)  gives 


5 (r  v.)‘ 

-3  h ♦ - V- 


, \2  6 

= or  (rlV]_)  - 


(D.21) 
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In  a close  neighborhood  of  r^  let 
6 5 (Ar)  1 

n 9 

u.  * c2  (Ar) 


(D.22) 


where  Ar  s r^-r.  Substituting  from  (D.22)  in  (D.20)  and  (D.21)  gives 
respectively 


-(nx+n2)  (a2-«4)  c^Ar 


n.+n_- 1 -vf  1 
12  o 


(Ar)"“l 


(D.23) 


nj+2n2-l 

”(n1<-2n9)  a3c1c22  (A r) 


= (arl)  - ‘ ~ (Ar)"1  (D.z!) 


ci (rivi ) 

"1 

Equating  indices  of  Ar  terms  and  solving  the  resulting  equations  gives 


n = 0,25 
1 


n2  = 0.50 


Combining  Equations  (D.22)  - (D.25)  gives 


6 


vf 


0.75(a4-or2 


j 1.25o'3r1 

) V (l-a1)v12 


1/2 


1/4 


(Ar) 


(D.25) 


u — - 
s 


jil-ai)  vx 
1 . 25o?3r1 


(Ar) 


1/2 


ujh/2  - a^6u 
U6  £ - (1^5)  6 


(D.26) 


The  numerical  algorithm  is  similar  to  the  one  described  in 
Appendix  B for  one -dimensional  model. 
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APPENDIX  E.  DYNAMIC  MODEL  OF  NSVR 


The  dominant  dynamics  of  a NSVR  element  are  expected  to  be  in 
the  vortex  chamber,  i.e.,  region  2.  The  inlet  and  the  central  orifices 
and  the  diffuser  are  expected  to  respond  much  faster  than  the  vortex  chamber. 
Hence,  steady-state  models  are  used  for  all  regions  except  region  2 of  NSVR. 
One-dimensional  approach  used  for  steady-state  modeling  is  also  used  for 
dynamic  model  development  as  described  below. 

Governing  equations  of  regions  1,  3,  and  4 of  NSVR  are  the  same 
as  those  given  in  Appendix  B.  However,  they  have  been  repeated  here  for 
easy  reference.  (See  Appendix  B for  assumptions  involved  in  modeling  of 
regions  1,  3,  and  4.) 


E.l. 

Region  1 

E.1.1. 

Governing  Equations 

Flow.  Q = CdlbKNy2(Pi-Pl)/p 

(E.l) 

Continuity:  u^=  -Q/2rTr^h 

(E.2) 

2 

Angular  Momentum:  pr^Q  2 

rr— — - pr.Qv  = 2nr1  hr  „ 

bhcdi  1 1 1 tl 

(E.3) 

where 

1.5  , 

Ttl  = 0.664m- vx  /VvXml 

(E.4) 

E.2. 

Region  2 

E.2.1. 

Assumptions 

below: 

Assumptions  in  addition  to  those  given  in  Appendix  B 

are  given 

1.  There  is  a time  delay  T.  and  a first  order  lag  T2  associated 
with  the  propagation  of  the  tangential  velocity  in  the  vortex  chamber. 

2.  Variations  of  through  flow  are  small  so  that  the  discharge 
coefficient  variations  are  negligible. 


E.2.2.  Governing  Equations 


Continuity: 


u 

2 


(E.5) 


Preceding  page  blank 
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Flow  Rate: 


2 CJ2nr2  [_2(P2"P3)  2] 

[2  3/2  1 

1 - |l-2 — ~ ;l  1 «... 

-(p2-p3)  + u22  J 


Region  4 


E.4.1.  Governing  Equations 


X 2 

Pressure  Recovery:  PD  = P3  + ^ P u3  + v3^ 


(E. 12) 


Equations  (E.l)  - (E.12)  can  be  linearized  (wherever  necessary) 
and  combined  to  result  in  a linear  dynamic  model  of  NSVR  which  can  be  used 
for  small  signal  AC  analysis.  Linearizing  equations  (E.l)  and  (E.2) 
gives  respectively: 


Ap1  - Lp1  = al  AQ 


Au^  = a^Q 


(E.13) 


(E.14) 


where 


al  = V2P(Pi*  - P!*)/Cdlbh 
a2  = -l/2rrr1h 

combining  Equation  (E.3)  and  (E.4)  and  linearizing 
Av^  * a3AQ 


(E. 15) 


where 


” 2pr  j( 

S 


Prlvl’ 


2 ^ 1 

1.992nr1zh  u . + pr^Q* 

L 1 ' v*ml  1 


Linearizing  Equations  (E.5)  and  (E.6)  gives  respectively 
Au2  * a^  Au3 


(E.16) 


♦Asterisk  refers  to  values  of  the  variables  at  the  operating  point. 


101 


Av2  = a 


Avl  + a6AUl 


5 1+T2D 

«<here 


--  


sign  |l -v2*2/C2J  (v2*/C2j  J 

2(p2*-p3*)  0 

C2S p +u2*‘ 

C3  = Cd2  Trr22/v2*2 

C4  = | 1 - v2*2/C2  | . 

Linearizing  Equation  (E.12)  gives: 

Ap3  " APo  = a12Au3  + a13^v3  (E .22) 

where 

a12  = - pug*/2 
a13=  - pv3*/2. 

Equations  (E.7),  (E.13)  - (E.22)  can  be  combined  to  give 


&P±  - AP0 
AQ 


_ T,D 

e,  + e2  e 1 

i+t2d 


(E.23) 


where 

el  = ai  +311  + a8a13  + a2  ^k2  + a4a9  + a6a7a13  + a6al0  " a4a12^ 
e2  = a3  ikl  + (a7a13  + aio)  ®5  } ' 

E.5.  Loading  Resistor 

For  an  orifice  type  of  loading  resistor  the  flow  rate  is  given  by: 
Q = Cd4A  \l2  pQ/p  (E.24) 

assuming  the  downstream  end  of  the  orifice  is  at  atomospheric  pressure. 
Linearizing  Equation  (E.24)  gives  (assuming  AQ  is  small  enough  so  that 
Cd3  remains  constant): 
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n . 


